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O For an arbitrary generalized quantum integrable spin chain we introduce a "mas- 

ter T-operator" which represents a generating function for commuting quantum 
transfer matrices constructed by means of the fusion procedure in the auxiliary 
space. We show that the functional relations for the transfer matrices are equiva- 
lent to an infinite set of model-independent bilinear equations of the Hirota form 
t-h for the master T-operator, which allows one to identify it with r-function of an 

integrable hierarchy of classical soliton equations. In this paper we consider spin 
chains with rational GL(A r )-invariant i?-matrices but the result is independent of a 
^ particular functional form of the transfer matrices and directly applies to quantum 

integrable models with more general (trigonometric and elliptic) i?-matrices and to 
supersymmetric spin chains. 
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1 Introduction 

The aim of this paper is to make precise the long anticipated connection between transfer 
matrices of quantum integrable models and r-functions of classical integrable hierarchies 
of non-linear partial differential equations. We show that properly defined generating 
function for commuting Hamiltonians of quantum models is a r-function in the sense 
that it obeys exactly the same bilinear identities of the Hirota form [H [2] as the classical 
r-function does. Since the operator-valued generating functions commute for all values 
of the auxiliary parameters (one of which being the familiar spectral parameter), there 
is no problem with their ordering in non-linear equations. 

A similarity between quantum transfer matrices and classical r-functions was first 
pointed out in [5] (see also [4]), where a discrete integrable dynamics in the space of 
integrals of motion of a quantum integrable model was introduced. This classical dy- 
namics was identified with the discrete 3-term Hirota equation with special boundary 
conditions. It is equivalent to the so-called T-system which is a system of fusion rela- 
tions among mutually commuting transfer matrices for bosonic [5] or supersymmetric [6] 
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quantum integrable spin chains. Solutions to the Hirota equation carry full information 
about spectral properties of the quantum model. The diagonalization of transfer matri- 
ces by means of the nested Bethe ansatz technique was shown to be equivalent to an 
"undressing" chain of Backlund transformations for the discrete Hirota equation. Later 
this approach was extended to supersymmetric integrable models [7]. 

An important further step was recently made in [8] , where an operator realization of 
the Backlund flow describing the "undressing" process was suggested and a commuting 
family of transfer matrices (T-operators on different levels of the nested Bethe ansatz) 
depending on a number of discrete variables (labels of representation in the auxiliary 
space) and on the spectral parameter was constructed. These T-operators obey the 
discrete Hirota equation not only for usual but also for supersymmetric integrable models, 
as it was directly demonstrated in [9] for GL(M|iV)-invariant spin chains. 

The question which remained open after these works was how to embed these T- 
operators depending on discrete variables into an infinite classical integrable hierarchy 
with continuous time flows compatible with the discrete ones. In the present paper we 
suggest an answer which is rather simple and natural. Below in the introduction we 
describe it in a concise form. 

In this paper we consider generalized quantum integrable spin chains with rational 
GL(iV)-invariant i?-matrix [10] 

N 

R x (u) = ul®l+ 6 Ji ® ( L1 ) 

which acts in the tensor product of the iV-dimensional vector representation 7Tn and an 
arbitrary finite dimensional irreducible representation tt\ of the group GL(N) labeled by 
a Young diagram A. Here are standard generators of the algebra gl(N), I is the unit 
matrix and u G C is the spectral parameter. A family of commuting operators (quantum 
transfer matrices or simply T-operators) , shown in figure [TJ can be constructed as 

T\u) = tr WA (R\u - &) <g> . . . ® R\u - £i)(I 0L ® 7r A (<?))) , (1.2) 

where are arbitrary parameters (inhomogeneities at the lattice sites) and g G GL(N) 
is called the twist matrix. The tensor product is taken over the quantum space (in the 



first space of (1.1 )). The trace is taken in the auxiliary space where the representation n\ 
is realized. The T-operators act in the physical Hilbert space of the model H = (C N )® L 
which is the tensor product of L local spaces C N of the vector representations. Formally, 
our setting includes also models with higher representations at the sites because they 
can be obtained by fusing several vector representations at the sites with properly chosen 
parameters The Yang-Baxter equation for the i?-matrix implies that the T-operators 
with the same g commute for all u and A and can be diagonalized simultaneously. They 
obey some functional relations which are given by the Cherednik-Bazhanov-Reshetikhin 
(CBR) determinant formulas [TT| [12] : 

T\u) = (JJ T\u-k)\ det ; T^- i+j \u-j+l), (1.3) 



k=i 



where A' x is the height of the first column of the Young diagram A and is the empty dia- 
gram. These relations mean that the transfer matrices for arbitrary diagrams functionally 
depend on the transfer matrices corresponding to the 1-row diagrams. 
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Figure 1: T-operator for an integrable, inhomogeneous spin chain with twisted boundary 
condition. 



We call such models "spin chains" in a rather broad sense, not implying existence of 
any local Hamiltonian of the Heisenberg type (integrable local interactions in general do 
not exist for inhomogeneous spin chains). However, even in the general case of arbitrary 
parameters & the model still makes sense as a generalized spin chain with non-local 
interactions. The "spin variables" are vectors from the spaces at each site. Instead, 
one may prefer to keep in mind integrable lattice models of statistical mechanics rather 
than spin chains as such. In either case the final goal of the theory is the diagonalization 



of the transfer matrices of the type (1.2) which is usually achieved by the nested Bethe 
ansatz method in one or another form (see, e.g., |13j). 

The key point of our approach is dealing with all the T-operators simultaneously by 
introducing their generating function of a special form which we will call the master T- 
operator. This generating function depends on an infinite number of auxiliary parameters 
t = {ti, t 2} ■ ■ ■} which we call time^j 

T(«,t) = ^ SA (t)T A ( tt ). (1.4) 



Here s\(t) are the standard Schur functions (2.13) and the sum goes over all Young dia- 



grams A including the empty diagram. By construction, the master T-operators commute 
for different values of the times: 

[T(M),TK,t')] = 0. (1.5) 

Our claim is that the master T-operator is the r -function of the KP hierarchy with 
the times £1,^2, • • • • In other words, the transfer matrices T (u) are Plucker coordinates 
for the element of the infinite dimensional Grassmann manifold corresponding to the 
r- function T(u,t). This means that T(u,t) obeys the Hirota bilinear equation 

(z 2 - z 3 )T (u,t+ [z^ 1 ]) T (u,t+ [z^ 1 ] + [z^ 1 )) 
+ (z 3 - Zl )T (u,t+ [z 2 1 ]) T (u,t + [zf 1 ] + *]) (1-6) 
+ ( Zl - z 2 )T (u, t + [Z3- 1 ]) T(u,t+ [^r 1 ] + fa 1 ]) = 



1 In a different form such an operator was introduced in [8] 
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with respect to the t- variables with fixed u. Here the standard notation t + [z -1 ] = 
{ti + z -1 , t 2 + \ z~ 2 , £ 3 + | z~ 3 , . . . } is used. This functional relation encodes infinitely 
many partial differential equations for T obtained after expanding the left hand side in 
inverse powers of Zj. 

Moreover, when one incorporates the it-dependence of the master T-operator, the KP 
T-junction extends to the r -function of the modified KP (MKP) hierarchy with the times 
to = u, ti, ti, This means that T(u, t) obeys another Hirota bilinear equation 

z 2 T (u + 1, t + [zr 1 ]) T(u,t+ \z 2 1 }) - ziT (u + 1, t + [z 2 1 ]) T(u,t+ [zf 1 ]) 

(1.7) 

+ ( Zl - z 2 )T (u + l,t + [z- 1 ] + [z^ 1 ]) T(u,t) = 

which encodes infinitely many differential-difference equations of the MKP hierarchy. In 
fact this statement is very general and is independent of a specific functional form of 
T x {u). It also does not depend on particular features of the quantum models in question: 
they can be lattice or continuous, bosonic or supersymmetric, with rational, trigonometric 
or elliptic /^-matrices, etc. 

Most part of the paper is devoted to the proof of this statement. Our technical tools 
are the realization of the KP hierarchy as a dynamical system on an infinite dimensional 
Grassmann manifold [H] and the free fermionic construction of the r-function [To*l ITB] 
developed by the Kyoto school. The main idea of the proof is to identify the generating 



series (1.4) for the master T-operator with the Schur function expansion of the r-function 



r n (t) of the MKP hierarchy: 

r n (t) = ^s A (t)c A (7i). (1.8) 

A 

The key fact proven in section 3 with the use of the free fermionic formalism is that the 
coefficients c\(n) (Plucker coordinates) obey the determinant relations of the form 

a;-i 

c aW = ( TT c ®( n - k )) det C( At _ i y)(n-j+l). (1.9) 
\.. ' i.7=l,...,A; 



k=l 



They mean that the Plucker coordinates for arbitrary diagrams are expressed through 
the basic ones, corresponding to the 1-row diagrams. Identifying n with the spectral 



parameter u and c\{n) with the T x (u), one can see that (1.3) coincides with (1.9). So, 
the statement that the master T-operator is the r-function appears to be equivalent to 
the Cherednik-Bazhanov-Reshetikhin (CBR) determinant relations for T x (u). 

For spin chains with rational i?-matrices, one can say more. In this case, the master 
T-operator has an explicit realization in terms of the co-derivative introduced in [9]. It 
was also remarked in [9] that the basic bilinear identity for the g/(M|iV)-character^]and 
their co-derivatives, which was the key part of the direct proof of the CBR formulas, 
might be related to the Hirota equation for a classical integrable hierarchy. Indeed, now 
we can see that the "master identity" from [HJ is a particular case of the bilinear equation 



(1.7) for the master T-operator. Therefore, the master T-operator contains all Baxter 



2 known as a particular solution to the discrete version of the KdV equation 
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Q-operators and T-operators on all levels of the nested Bethe ansatz. Presumably, this is 
also true for models with other types of i?-matrices. The most important characteristic 
property of the master T-operators for finite spin chains with rational i?-matrices is that 
they are polynomials in u: 

L 

T(u,t)oc JJ( u — Uj(t)). (1-10) 

The class of r-functions with this property is well studied in the theory of soliton equa- 
tions and can be characterized [J71 HE] in terms of algebraic geometry of curves with 
cusp singularities. Furthermore, the dynamics of the roots Uj as functions of times £j, is 
known [19] to be given by the rational Ruijsenaars-Schneider model. 

Sections 2-4 of the paper contain the necessary material on the classical KP and MKP 
hierarchies and their r-functions. Although the main result can be proved without using 
the formalism of free fermions, we found it instructive to present all arguments in the 
fermionic operator language of the Kyoto school. In section 2 the free fermion operators 
and the vacuum states are introduced and the vacuum expectation values are defined. 
In section 3 the r-function is introduced as the vacuum expectation value of a group-like 
element and bilinear identities for it are derived. In most cases we use the standard 
notation and conventions from [Tol ITS] . Section 4 is devoted to a detailed review of 
rational solutions to the MKP hierarchy. It is precisely the class of solutions relevant to 
quantum integrable models with rational i?-matrices. The core of the paper is section 
5, where the master T-operator is introduced and shown to satisfy the classical bilinear 
identities and Hirota equations. 

2 Free fermions 

In this section we recall the formalism of free fermions [T51 (TB] used for construction of 
r-functions of integrable hierarchies. 

2.1 Fermionic operators 

Let ■*/>„,■?/>*, n e Z, be free fermionic operators with usual anticommutation relations 
[i> n ,i>m]+ = VPn^*m\+ = °> H'n,^m}+ = ^mn- They generate an infinite dimensional 
Clifford algebra. We also use their generating series 

fcez fcez 

which can be regarded as free fermionic fields in the complex plane of the variable z. 
Of particular importance is the operator 

J+ = ^t fc Jfc, Jh = ^2Mj+k> ( 2 - 2 ) 

fe>i jez 
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where tk are arbitrary parameters (called times). It is convenient to denote the set of 
times by t = {ti, t 2 , ■ ■ ■} and to introduce their generating function 

£(t, = 

k>l 

It is easy to check that the fields ip(z) : ip*(z) transform diagonally under the adjoint 
action of the operator e J + : 

e J+ iP{z)e- J+ = e f(t ' 2 V(^) , e J+ *fj*{z)e- J+ = e~^ z ^\z). (2.3) 

In terms of the symmetric Schur polynomials h^{t) defined by 

e «M) = Y^h k {t)z k (2.4) 

k>0 

the corresponding formulas for ip n ,ip* can be written as 

e J +^ n e- J + = J2^n- k h k (t) , e J +<e- J + = ^< +fc /i fe (-t). (2.5) 

fc>0 fc>0 

2.2 Dirac vacua and excited states 

Next, we introduce a vacuum state |0) which is a "Dirac sea" where all negative mode 
states are empty and all positive ones are occupied: 

^„|0) = 0, n<0; C|0) = 0, n > 0. 

(For brevity, we call the modes with indices n > as positive). With respect to this 
vacuum, the operators ip n with n < and ip* n with n > are annihilation operators while 
the operators ip* with n < and ip n with n > are creation operators. The dual vacuum 
state has the properties 

(0|< = 0, n<0; (0|Vn = 0, n > 0. 

We also need the "shifted" Dirac vacua \n) and the dual vacua (n| defined as 

{■0„_i . . . ^1-00 |0> , n>0 
C---C 2 Ci|0), n<0 

r (oi^svi ---C-i, «>o 

<n| = { 

[ (0| ^-1-0-2 •••'0n, W < 

Similarly to the properties of the Dirac vacuum, for the shifted Dirac vacuum we have: 

ipm \n) = 0, m < n; ip* m \n) = 0, m > n, (2.6) 
(n| ip m = 0, m > n; (n\ ip* m = 0, m < n. (2.7) 
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We will also use 



tP n \n) = \n + l), <|n+l) = |n>, (2.8) 
(n + l\if) n =(n\ (n\il>* = (n + l\. (2.9) 

Excited states are obtained by filling some empty states and creating some holes. Let 
us introduce a convenient basis of states in the fermionic Fock space J 7 . The basis states 
|A,n) are parameterized by Young diagrams A and numbers n of the Dirac vacua in the 
following way. Given a Young diagram A = (Ai, . . . , A^) with I = £(X) nonzero rows, 
let (a|/5) = («i, . . . , (Xd(X)\Pi, ■ ■ j 0d{X)) be the Frobenius notation for the diagram A [20J. 
Here d(X) is the number of boxes in the main diagonal and = Aj — i, = A^ — i, where 
A' is the transposed (reflected about the main diagonal) diagram A. In other words, a« 
is the length of the part of the z-th row to the right from the main diagonal and is the 
length of the part of the z-th column under the main diagonal (not counting the diagonal 
box). Then 

| A, n) := VC-A-l • • • C-/3 d(A) -l ^n+a d(x) ■ ■ ■ *Pn+ ai \n) 

(2.10) 

(A, n\ := (n\ ip* n+ctl . . . C +Qd(A) ^n-p d{X) -i ■ ■ ■ VV-A-i- 

The definition of the vacua implies that e J+ \n) = \n) while the dual coherent states 
(n\ e J+ are expanded as linear combinations of the basis vectors as follows: 



e J+(t) = ^(-l) 6(A) ^(t)(A,n|, (2.11) 

A 

where the sum runs over all Young diagrams A including the empty diagram, 

rf(A) 

6(A) = 5>< + 1 ) ( 2 - 12 ) 

i=l 

and s\(t) is the Schur polynomial corresponding to the diagram A. It can be expressed 



through the functions hk(t) defined in (2.4) (which are Schur polynomials for one-row 



diagrams) with the help of the Jacobi-Trudi formula 

*A(t) = .. det V^(t). (2.13) 

2.3 The expectation values and normal ordering 

The vacuum expectation value (n| . . . \n) is a Hermitian linear form on the Clifford algebra 
defined on bilinear combinations of fermions by the properties (n\n) = 1, (n| ipiipj \n) = 
(n\ ?p*?p* \n) = for all i,j and 

(n\ ipiip* \n) = 5ij for j < n, (n\ i^iip* \n) = for j > n. 



One can see that the basis vectors (2.10) are orthogonal with respect to the scalar product 
induced by the expectation value: 

(A,n| /j,, m) = <W<W- 
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Expectation values of general products of fermionic operators are given by the Wick's 
theorem. Let Vi = ^jVijipj be linear combinations of ^-'s only and w* = Y^j w tj 1 i , *j be 
linear combinations of ip^s only, then the standard Wick's theorem states that 

(n\ vi . . . vmw* m . . . wl \n) = det (n\ ViW* \n) . 

i,j=X,...,M 

One may define the normal ordering *(...); with respect to the Dirac vacuum |0): 
all annihilation operators are moved to the right and all creation operators are moved to 
the left taking into account their anti-commutativity under mutual permutations. For 
example, \^x\ = ~Ml = tpi^i ~ 1, and '.tp^n'. = ^* m ^Pn - (0|V£A|0). More 
generally, for any linear combinations fo, fi, ■ ■ ■ , fm of the fermion operators t/>? we 
have the recursive formula 

m 

/o:/i/2 ...f m : = :/ /i/ 2 ..../;,: • » d j ' (o| ./;,/• |o) :a/ 2 . . . % . . . f m \ (2.14) 

where fj means that this factor should be omitted. 

In fact a normal ordering can be defined with respect to any vacuum state. Another 
useful normal ordering is the one defined with respect to the bare vacuum |oo), which is 
the absolutely empty state. We denote this normal ordering by * (. . .) * . It means that all 
■0*'s are moved to the left and all ?/>'s to the right, with taking into account the sign factors 
appeared each time one operator is permuted with another. For example, *?Pm4>nx = 
CV'm x^mx = which is equivalent to £Vv0m£ = ~ (oo| i) n ^ m |oo). In 



general, for this normal ordering equation (2.14) holds for the vacuum expectation taken 
with respect to the bare vacuum. 



2.4 Group-like elements 

Neutral bilinear combinations ^2 mn bmnip^ipn of the fermions, with certain conditions on 
the matrix b = (b mn ), generate an infinite-dimensional Lie algebra [16]. Exponentiat- 
ing these expressions, one obtains an infinite dimensional group (a version of GL(oo)). 
Elements of this group can be represented in the form 

G = exp(^M^fc) (2.15) 

j,feez 

The inverse element can be written in the same way with the matrix (— b mn ). An impor- 



tant example is provided by the operator e J+ (2.2). 

It is straightforward to see that the group elements of the form (2.15) obey a rather 
special property that the adjoint action of such elements preserves the linear space 
spanned by the fermion operators ip n and the dual space spanned by More pre- 
cisely, we have: 

i i 

where the matrix R = (R n i) of the induced linear transformation is given by R = e b . 
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Let us note, following the works of the Kyoto school [ISlET], that the group elements 
can be equivalently represented as normal ordered exponents of bilinear forms. For 
example, G given in ( |2.15[ ) can be written as 

G = * exp(^T B ik ijj*ipk) * = det(I+P + B) : exp(^ A ik ^k) '. (2.17) 



with the matrices A, B determined by the matrix b in (2.15) according to the formulas 

m 

B = e b -I, B-A = AP + B. (2.18) 

Here I is the unity matrix and P + is the projector on the positive mode space ((P + )ik = da- 
for i, k > and otherwise). The product of two group-like elements is also a group-like 
element: 



where C = A + B + AB. As a slight generalization of (2.17), one can obtain the formulas 

(n| x e E ifc s«^ ^ = det j 7 + B (n)j ^ ( 21Q ^ 

where B^ = (B ik ) 

i : k> n is the half-infinite matrix obtained from B by truncation up to 
the n-th row and n-th column and 

det (/ + (l + B^y 1 if r,s>n 

(n|^*^e EifcBifc ^ fc ^ r \n) = { " r ' s (2.20) 

otherwise. 



Note that the r.h.s of the last formula is the (rs)-minor of the matrix / + B^ (up to a 
sign) . 

The general theory of normal ordering and the group elements in the the Clifford 
algebra was given in [21]. What is important here, is that the normal ordering allows one 
to represent in the form (2.17) not only the group elements but also some non-invertible 
elements of the Clifford algebra that satisfy the same commutation relations (2.16) with 
some matrix R. 

We call the elements G of the Clifford algebra, such that the commutation relations 



(2.16) hold with some (not necessarily invertible) matrix R, the group-like elements. If 
the matrix R fails to be invertible, so does G, as an element of the Clifford algebra. In 
this case it can not be represented in the exponential form (2.15). However, in general 
it still admits a representation as a normal ordered exponent of bilinear forms in the 
fermionic operators. 

Let us give an example. Let <3>* be arbitrary linear combinations of the fermion 
operators ip n , ip*, respectively, and consider the element 

G = e^** = * e a *** * = 1 + = 1 + a 7 - a W, 

where 7 := (00 1 |oo) and a, (3 are related as e 7/3 = 1 + 7a. For general values of 
a (and for 7 7^ 0) this element is invertible and the two representations are equivalent. 
However, at a = — 1/7 the invertibility breaks down and the element G becomes 

\T/ <T>* 

_ 

(oo|^$*|oo)' 
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which can not be written in the exponential form (2.15) but can be represented as the 
normally ordered exponent. 



From (2.16) it easily follows that any group-like element obeys the commutation 
relation 



(2.21) 



fcez fcez 
which is equivalent to the bilinear identity 

E w ^ G \ y ) w\ ^ G \ v ') = E < u \ G ^ \ y ) (^i G ^ \ v ') 



(2.22) 



valid for any states \V) , \V), (U\ , (U'\ from the space T and its dual. This is the basic 
identity for deriving various bilinear equations for the r- function (see below). 



2.5 The generalized Wick's theorem 



Let Vi = J2j v iji J j be linear combinations of ^-'s only and w* = ^2jW*jip* be linear 
combinations of ?/>*'s only, as before. A more general version of the Wick's theorem given 
in [15] is 



(n 


Vi... V m 




...w{G 


n) 


(n 


G 


n) 



det 

i,j=l,...,m 



(n 




w*G 


n) 


(n 


G 


n) 



(2.23) 



where G is any group-like element. Here we imply that {n\ G \n) ^ for all n but in fact 
this restriction can be removed by multiplying the both sides by an appropriate power of 
(n| G \n). The similar formula with the exchange v j <-> Wj also holds. This statement can 



be proved by induction using the bilinear identity (2.22) (see Appendix A). The following 
particular case is often useful in applications. Writing (n — m\ = (n \ ip n -i ■ ■ ■ ipr, 



m > 0, we get from (2.23) 



(n — m 


u 


* 

m " 


..w\G 


n) 


(n 


G 


n) 



det 

i,j=l,.... 



(n 


Ipn- 




*G 


n) 


(n 


G 


n) 



with 



(2.24) 



For our purpose, it is important to note that there exists an alternative determinant 



representation of the expectation value in the l.h.s. of (2.24), which is another form of 
the generalized Wick's theorem: 



(n — m 


u 


* 

m ' 


..w\G 


n) 


(n 


G 


n) 



det 

i,j=l,.... 



(n-j 


w*G 


n-j+1) 


(n-j+l\G 


n-j+1) 



(2.25) 



Indeed, the first columns of the two matrices in the r.h.s. of (2.24) and (2.25) coincide 



and one can show that the j-th column of the matrix in (2.24) is equal to the j-th column 



of the matrix in (2.25) plus a linear combination of the first j — 1 columns (see Appendix 



A). We also note an equivalent form of (2.25): 



(n 


Gv x .. 




n — m) 


(n 


G 


n 


) 



det 

i,j=l,...,m (n 



n- 



+ Gvj \n-j) 
-j+l\G\n-j+l)' 



(2.26) 



which is obtained from it by conjugation of the states and operators. 



12 



3 The r-function and the Baker- Akhiezer functions 



As it has been established in the works of the Kyoto school, the expectation values of 
group-like elements are tau-functions of integrable hierarchies of nonlinear differential 
equations. In particular, 

r (t) = (0| e J +^G |0) (3.1) 

is the tau-function of the KP hierarchy depending on the times t = {£1,^2, • • •}, which 
means that it obeys an infinite set of Hirota bilinear equations [H [2] . More generally, 
introducing a "discrete time" n via 

r„(t) = (n| e J+{t) G \n) (3.2) 

one obtains a r-function which solves the modified KP (MKP) hierarchy [221 E31 12/ 1 '' 



Equations of the MKP hierarchy are differential-difference equations which include shifts 



of the variable n. At each fixed n, the r-function (3.2) is a r-function of the KP hierarchy, 



so the n-evolution represents an infinite chain of Backlund transformations. The meaning 
of the discrete variable n depends on the class of solutions. For some classes of solutions 
n can be regarded as a continuous, or even complex, variable. The last point will be of 
a prime interest for us since in section [5] n will be identified with the spectral parameter 
u of the quantum transfer matrix of a spin chain. 

3.1 Bilinear equations for r-function 
3.1.1 The bilinear identity 



The r-function obeys a bilinear identity which is a direct consequence of (2.22). Setting 



\V) = |n), \V) = \n') with n > n' , where \n) and \n') are two shifted Dirac vacua, we 
have 

J2(U\^kG\n) (U'\^* k G\n'} = (3.3) 



because either ipk or ipl kills the state in each term in the r.h.s. of (2.22). Now, setting 
(JJ\ = (n + 1| e J+(t) , (U'\ = (n' - 1| e J+{t '\ we can write 

= ^(n + l|e J+(t ^ fc G|n)(n / -l|e J+ ( t VfcG ! |n / ) 



res. 



z- 1 (n + 1| e J+{t) i){z)G \n) (ri - 1| e J + ( - t ' ) ip*{z)G \n') 



res. 



-t',z) z n-n' ^| e J + (t-[,-i]) G | n ) ^/| e J + (t'+[,-i]) G | n /j 



where we have used the commutation relations of the fermion operators with e J+ ^ and 
the bosonization rules 

(n\ ^)(z)e J+{ - t) = z n ~ l (n - 1| e^^-l*" 1 ]) 

(3.4) 

(n\ i)*(z)e J+{t) = z~ n (n + 1| e J+(t+[2 1]) 



3 In [21] a slightly more general version, called there a coupled modified KP (cMKP) hierarchy, was 
considered. 
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established by the Kyoto school. Here and below we use the standard notation 



t ± [z- 1 ] = {ti ± z-\ t 2 ± ± ■ ■ •}• 

In this way we arrive at the bilinear identity 

j£ e «t-t',*)s»-»V n (t - [z-^Tnit + [z- l })dz = (3.5) 

for the r-function valid for all t,t' and n > n' . It encodes all the PDE's of the MKP 
hierarchy. 

The choice of the integration contour C depends on the type of the time evolution. 
Formally, the evolution factor e^ t_t ' z )z n ~ n has an essential singularity at oo, and the 
contour should then encircle the whole complex plane (only the singularity at oo stays 
outside the contour). This is the usual prescription and if only a finite number of times 
tfc are nonzero, then it is correct. However, in general, when the values of the times 
are such that the factor e^ t_t ' ,2 ^2; n_n ' has singularities for finite values of z, then the 
prescription is as follows: the contour C must encircle all singularities of the function 
r ra (t — [^ _1 ])r n /(t / + [z -1 ]), leaving all the singularities of e^ t_t ' z )z n ~ n outside the contour. 

Remark. There is a freedom to multiply the r-function of the MKP hierarchy by 
an exponent of any linear form of times with constant coefficients and by an arbitrary 
function of n: 

r n (t) -> C(n) exp Cfct fc ) r n (t). (3.6) 

k>l 

Clearly, this transformation preserves the form of the bilinear identity. By noting that 



exp 



^djipjipj * \n) = C(n) \n) , d n 



C{n) 



C(n + 1) 



we see that the transformation r n (t) — > C(n)r n (t) means G ->■ G*eZi d ^* for the 
group-like element. In the MKP hierarchy, there are no restrictions on the form of 
C(n). Such restrictions appear in the more general 2D Toda lattice hierarchy, where this 
function is constrained to be of the form C(n) = Ce c ° n with constant C, Co- 

3.1.2 3-term bilinear equations 



Setting n' = n and t' k = tk — \{ z \ k + z 2 k + z z k ) m (3-5), we see that the essential 
singularity at oo splits into 3 simple poles at zi,z 2 ,z 3 : 



e £(t-t',z) 



'1 - ^)(1 - ^)(1 - ^) 



Taking the residues, we arrive at the 3-term relation 

(z 2 - z 3 )r n (t - K 1 ]) r n (t - [z^\ - [z^ 1 ]) + (231) + (312) = 0, (3.7) 
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where the last two terms are obtained from the first one by the cyclic permutations of 
the indices. Another choice is n' = n and t' k = t k + \z^ k — \\z7 + z7 k + z7 k ), then 



3 «(t-f,*) 



u k — "fc -r k « k \^i 
1 - ^ 



This leads to a slightly more general equation 
(z - Zl )(z 2 - z 3 )r n (t - [z l ] - [z7 1 ]) r n (t - [z^ 1 ] - [z^ 1 ]) + (231) + (312) = 0. (3.8) 

Setting n' = n — 1 , t' k = t k — \{z7 + z 2 k ), we see that the essential singularity at oo 
splits into simple poles at z\, z 2 : 



ze 



«t-t',*) 



Besides the residues at these points, there is also a contribution from the residue at oo, 
so we obtain the 3-term relation 



Z 2 T n+1 (t - [Z 2 X ]) T n (t - [z x X ]) - Z X T n+X (t - [Z x X ]) T n (t - [z 2 *]) 



+ (Zl ~ Z 2 )T n+ i(t)T n (t - [z x l ] - [Z 2 *] 



(3.9) 



0. 



It can be formally regarded as a particular case of (3.8) in the limit z 3 — > oo, z — > 
0. The limit z 3 — > oo is smooth while the other one requires to assign a meaning to 
lim r„(t — [zq ]). The form of the evolution multiplier, z n e^' z \ sug gests to set formally 

2o->0 



lim r n (t ± [z l ]) = r„ T i(t) 



(3.10) 



which actually holds if the function r n (t — [z~ 1 ]) for arbitrary n can be analytically 
continued from a neighborhood of oo to a neighborhood of the point z = and is regular 
there. If it is singular, as in the case of solutions relevant to quantum spin chains, (3.10) 
should be substituted by a more general prescription (for an example see ( 3.14[ ) below) 
which, however, also allows one to regard (3.9) as a particular case of (3.8). 



3.2 Examples of r-functions 

Here we consider some particular cases of the general fermionic construction of r-functions 
which will be important for us in section 5. 

Let us start from a very simple example. Set \1/ = ip(p) + oti}j{q) and consider the 
operator •e l3q "^*^ • (with |r| < min (\p\, \q\)). It is a group-like element for any (3 but at 
/? = A) = ((0| $fip*(r) |0)) _1 it is not invertible and equals (3 ^ip*(r). Therefore, 

r n (t) = (n|e J+(t) ^*(r) \n) 

is a r-function of the MKP hierarchy. Then we would like to consider the limit r — > 0. 
The limit is singular and requires some care. From 

r n -V*(r)|n) = ^r n -'~V/>) = ^"'"Vri") = |n-l) + G(r) 

i i<n-l 
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we see that in order to get a well-defined limit, one should multiply the r-function by 
r n-i De f ore tending r — >■ (this is just a transformation of the form (3.6)). Then we 
come to the conclusion that 



(t) = (n\ e J+(t) ^ \n - 1) = (n| e J+(t) U(p)+aip{q)) \n-l) = p^e^ + aq^e^ 



is a r-function. Indeed, it is the r-function for the 1-soliton solution of the MKP hierarchy. 

This example can be generalized. Fix a finite set of distinct points pi G C, % G /, and 
construct a (finite) linear combination of the operators ip(z) and their ^-derivatives at 
the points pf. 



Examples of the operators are ?p(p), ipip) + aip'(p), ?P(pi) + cnA(P2), ip(pi) +otip'{j>2) + 
P4 ! "{P2) etc, where ip'(z) = d z ip(z). Clearly, (n\ e J+ ^> 'e^ 1 ^*^ m m \n) is a r-function for 
any fixed (3. Let us take (3 = ((0| tyjip*(r) |0)) _1 , then the group-like element is non- 
invertible but (n\ e J+<yt ^ iip*(r) \n) is still a r-function. The limit r — > can be performed 
as in the previous example and we conclude that 



r n (t) = (n| e J+{t) ^/ \n - 1) 



is a r-function. 



This example can be further generalized by considering products of operators of the 
form Take N operators of the form ^/i for different sets Ii, . . . , In- ^h, • • • , ^i N - As 
before, we can construct the group-like elements 'e 13 ^ 1 ^ and consider a r-function 
(n\ e J+(t). e PN*i N r(r N ). . . . .gftSx^Cn). Note that if the auxiliary points r, are all 
distinct from the points labeled by the sets Ij, rj ^ p^ for all j, k G /j U . . . U I^r, 
then these operators commute. Again, we choose {3j = ((0| ^i.i/;*(rj) |0)) _1 , then each 
operator becomes the product Pj^i.i/j*(rj) (non-invertible). Thus the above r-function 
reduces to a r-function of the form ((3i ■ ■ • Pn) (n\ e^^^/j • ■ ■ ^i N tp*(rN) ■ • •'0*( r i) \n). 
In order to implement the limit rj — > 0, we redefine rj — > erj with e — > 0, then it is not 
difficult to verify that 

£ (n-i)jv-|iv(iv-i)^*( £rjv ) ^*(en) |n) = (ri.-.r^-^Ajv^) \n-N) +0(e) 
e -nN-\N(N-i)^ erN ^ ^(en)|n> = (n . . .r N ) n A N (n) \n + N) + 0(e), 

(3.12) 

where 

A N ( ri ) = det r 3 ' 1 = TT(ri - j-,-) 

i,j=l,...,N -*- -*- 

%>j 

is the Vandermonde determinant. The first equation says that in order to get a well- 
defined limit, one should multiply the r-function by £ ( n - 1 ) N ~2 N ( N ~ lS > before tending e — > 



0. The r r dependent factors (ri...rjv) n+1 A7v(rj) in the r.h.s. of (3.12) are irrelevant 



because they can be eliminated by a transformation of the form (3.6). We thus conclude 
that 

r n (t) = (n\ e J +^ 7l . . . ttjjn - N) (3.13) 

is a r-function. 
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In this case the prescription (3.10) is not directly applicable. The behavior of the 



function r(t ± [z *]) at Zq — > can be found by using the bosonization formulae (3.4) 
in the opposite direction and moving the fermion operators ip(zo) or ip*(zo) to the very 



right position. In this way we obtain a more general version of the prescription (3.10): 



lim [(-zo^T^tTizo 1 ))} =r n ±i(t) 

2o->0 



(3.14) 



which is equally enough to deduce the bilinear equation (3.9) from (3.8) as a limiting 
case. 



3.3 Pliicker coordinates 

One may expand the r-function in the Schur polynomials: 

T n(t) = ^c A (n)s A (t). 



(3.15) 



The coefficients c A (n) ("Pliicker coordinates") can be determined from the formula (3.2) 
by inserting the complete set of states in between the operators e J+ and G and using 



(2.11): 



so 



r n {t) = H eJ+W l A > n) (A, n\ G \n) = ^(-l) fc(A) s A (t) (A, n\ G \n) 

A A 

c A ( n ) = {-l) b ^ (X,n\G\n) 

= (-l)KA) ( n | ^ . . . ^^Vv-A^-i • . • ^n- h -iG \n) , 



(3.16) 



where 6(A) is defined in (2.12). The Schur function expansions of r-function are also 
discussed in 125, 261. 



The coefficients c\(n) generalize the characters of representations of the linear groups. 
Remarkably, they admit determinant representations of the Giambelli type as well as of 
the Jacob i-Trudi type which express c\{n) for arbitrary A through c\{n) for the hook and 
the one-row diagrams, respectively. 



3.3.1 Formulas of the Giambelli type 



Applying to (3.16) the Wick's theorem in the form (2.23), we obtain the Giambelli-like 



formula for C\(n) 



c A (n) = (-lf x \(n\G\n})- d ^ +1 ^ det M C^V^-iG |n> 

z,j=l,...,d(\) 



(c (n))- d W +1 det ^(n). 

i,j=l,...,d(A) 



(3.17) 



Here C( Q , i | i g.)(n) are the Pliicker coordinates corresponding to the hook diagrams (oti\f3j 
(ctj + 1, 1^'). In the last expression, we have taken into account that 



(n\ G \n) = r n (0) = c (n). 
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Formulas of the Giambelli type for the expansion coefficients, or Pliicker coordinates of 
the r-function were given in [25] • Using the Jacobi identity for determinants, it is easy 
to see that they are equivalent to the 3-term bilinear relation for the coefficients c\ given 
in M\. 



3.3.2 Formulas of the Jacobi-Trudi type 



Alternatively, we may apply to (3.16) the Wick's theorem in the form (2.25). The details 



are given in Appendix B. The result is 



c\{n) 



(3.18) 



where c s (n) := C( s -i\o)( n ) — ( n ~ 1| VCfs-i^ \ n ) are the expansion coefficients for one-row 
diagrams and X[ is the height of the first column of the diagram A. 

Therefore, the Pliicker coordinates of any MKP r-function satisfy the determinant 



relations (3.18) of the Jacobi-Trudi type sometimes called quantum Jacobi-Trudi for- 



mulas. Note that the transformation c\(n) — > C{n)c\(n) with arbitrary C(n) preserves 



the determinant formula (3.18). Clearly, this freedom corresponds to the possibility of 



multiplying the r-function by any C(n), see (3.6). 



In fact the inverse statement is also true: given any set of quantities c\(n) that satisfy 
the determinant relations (3.18), the function c\(n)s\(t) is a r-function of the MKP 



hierarchy. In order to prove this, it is enough to show that for any given set of quantities 
co(n) = C0(n), c s {n) = C( s _i| )(n), s > 1, one can find a (n-independent) group-like 
element G such that 



,(n + 1) = (n| ipn +s G \n + 1) = (n\ il)* n+s G%l) n \n) , n G Z, 



0,1,2, 



(3.19) 



Let us give a sketch of proof. First of all we note that the element G, if it exists, is not 



unique because * e^> 



n) = \n) for any strictly lower triangular matrix A. This 



means that we can consider only upper triangular matrices. For simplicity we assume 
that c$(n) 7^ for all n. Then we can (and from now on will) put Co(n) = 1 without 
any loss of generality because this is the matter of normalization due to the freedom to 
multiply the MKP r-function by an arbitrary function of n as in (3.6). It is convenient 
to formally put Ck{n) = for k < 0. Consider the infinite matrix C nm = c m _ n (n + 1) 
(m,n G Z). It is an upper triangular matrix with l's on the main diagonal. One can see 
that for any matrix C of this form there exists a strictly upper triangular matrix B such 
that 

C nm = (n\ * exp(j^ B ik ^kj * |n + 1) . 



i<k 



Indeed, according to (|Z20[, C nm = [(I + B^y 1 ]^ 
strictly upper triangular B it holds M^ n ' = (I + B^) 
determined as B^ = (C _1 )ifc — 5ik- 



Set M = {I + B)-\ then for 
and so the matrix B is uniquely 
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3.3.3 Example: characters 



As an example, let us consider the r-function corresponding to the coherent states 



G \n) = exp I ^x k J- k J 
\fc>i / 



1 77.) 



where J-k are given by equation (2.2), Xk = \ tig k ) g G GL(N). The r-function (3.2) is 
explicitly given by 



N 



r n (t) = exp ( kt k x k ) = exp ( Yl tkP i I ( 3 ' 2 °) 

K k>l 



where Pi are eigenvalues of the matrix g. This r-function does not depend on n. There- 
fore, the coefficients c\(n) are n-independent and coincide with GL(iV)-characters eval- 
uated for the element g G GL(N): 



Xx{g) 



det 



s A(xj 



l<i,i<iV {pj 



Xi-i+N 



) 



n 



KKk<N 



(Pi - Pk) ' 



(3.21) 



The formulas (3.17) and (3.18) become the well-known determinant formulas for charac- 



ters [20] (Giambelli and Jacobi-Trudi respectively). 

Remark. This simple example may be somewhat misleading because in this case the 



r-function (3.20) is simultaneously the r-function with respect to the "negative" times 
t-k — Xk- For more complicated solutions constructed below from the group element 
g G GL(N) with Xk = rtr g k this is not true in general. 



3.4 Baker- Akhiezer functions and wave operators 

The Baker-Akhiezer (BA) function and its adjoint are given by the Sato formulas 



(n + 1 


e J +^{z)G 


n) 


(n 




n) 



Mt, z) = z ^ y^ » = v 17 (3.22) 



(n - 1 


e 


J +Mif>*(z)G 


n) 


*( 


n 


e J +^G 


n) 



«(M) = = ^z^_gjfj^v . n.m 



In terms of the BA function and its adjoint, the bilinear identity (3.5) acquires the form 



i/> n {t,z)i/)* n ,{tf,z)dz = 0. (3.24) 



c 



An important role in the theory |27j is played by the wave operator (or dressing 
operator) W(n, t) directly related to the BA functions. It is a pseudo-difference operator 
of the form 

W(n,t) = ^w fc (n,t) e - fca ", w (n,t) = l, (3.25) 



fc>0 
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and its formal inverse 

W-\n,t) = ^V w Si£(n+l,t), w* (n,t) = l. (3.26) 



fe>0 



The Baker- Akhiezer functions are obtained by applying the wave operators to the bare 
exponent z n e^ t,z ': 

V„(t, z) = W(n, t)z n e^ =(l + + + ..) z n e^ (3.27) 

\ z z 2 J 

V£(t,*) = (W- 1 (n-i,t))t«-» e -€( t '«) = ^ + + + . . }j z -n e -t{t,*) 

(3.28) 

(by definition, the adjoint operator for A = f{n)e k9n is A* = e~ k9n f(n)). Using formulas 



(3.22), (3.23) we obtain 



, . h k {-d)r n {t) h k {d)r n {t) 
w k {n,t) = — , w k {n,t) = 7- — , (3.29) 

T n (t) T n (t) 

where the standard notation d = {d tl , \d t21 |9 t3 , . . .} is used. These formulas allow us to 
express the wave operator and its inverse in terms of the r-function as follows (cf. [28J): 

W (n, t) = y ^(-^W e -*ft. (3.30) 

W-\n,t) = Ve- M " Mg^HlW (3.3i) 
z — ' r n+ i(t) 
fc>0 n+iV ' 

The (pseudo-difference) Lax operator of the MKP hierarchy is obtained by dressing 
the shift operator e dn with the wave operators: 

L = W{n, t)e 9 " JU -1 (n, t) = e dn + a {n, t) + ai {n, t)e~ dn + ... (3.32) 

The Lax operator obeys the Lax equations 

d tk L = [(L fc )> , L), (3.33) 

where (L fc )> means the (finite) part of the series containing the shift operators e^ dn with 
3 > 0. 



4 Rational solutions of the KP and MKP hierarchy 

In this section we study rational solutions of the KP and MKP hierarchies. For these 
solutions, the r-function is a polynomial function of the times n and (possibly multiplied 
by an exponential function of a linear combination of the times). Stressing the fact that 
n for solutions of this class is not necessarily integer but can be any complex number 
u e C, we change the notation n — > u. 
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4.1 The construction of rational solutions 



The general construction of rational solutions to the KP hierarchy was given by Krichever 
in [T7j, see also [18]. This construction can be directly extended to the MKP hierarchy. In 



this section we show how it works starting from the bilinear identity (3.5). The solutions 
are obtained in the Casorati determinant form which is a discrete analog of the Wronskian 
determinant (see, e.g. [29]). 

4.1.1 Baker- Akhiezer functions and bilinear identity 

First of all, let us specify the bilinear identity for the case at hand. Since in the solutions 
we are going to consider u is not necessarily integer, and oo are in general branch points 
of the BA function. In this case the appropriate form of the bilinear identity is 



^(t, 2 ;)^(t , ,z)rfz = 0, (4.1) 



C[0,c 



where [0, oo] is an arbitrary cut between and oo and the contour C[o i00 ] is such that it 
encircles the cut but does not enclose any other singularity of the BA functions. 

According to the Krichever's theory of rational and more general algebro-geometric 
solutions, they can be characterized and explicitly constructed by fixing certain analytic 
properties of the BA function on a Riemann surface. For rational solutions, the Riemann 
surface is the Riemann sphere (compactified complex plane), which represents a genus 
zero algebraic curve with singularities. Correspondingly, the BA function is, in this case, 
a rational function on the complex plane multiplied by power-like and exponential factors 



which give the required asymptotics (3.27) (the essential singularity at infinity). Let us 



assume the following ansatz for the BA function: 

Mt, ,) = ,V«W (l + ^ + . . . + 5M) , (4.2) 
where the coefficients Wi depend on the times u, tj. In fact this is just the truncated 



3 ' 

series (3.27). The multiplicity N of the pole at z = of the function z~ v e~^ t>z 'i(} u (t, z) 



is a discrete parameter characterizing the class of solutions to be constructed. 



Given the BA function of the form (4.2), an important dynamical information is 



contained in the adjoint BA function. In general, it may have poles of arbitrary order at 
some points pi G C (we assume that pi ^ 0). As we shall see soon, the bilinear identity 
is consistent when the number of these points does not exceed N, and, in case of general 
position their number is just equal to N. Let us adopt this and fix N distinct points 
Pi G C, pi 7^ 0, and also assume that the adjoint BA function has poles of orders Mi + 1 
at pi. Therefore, its general form is 

\ i=l m=0 V y% > J 

where a,i m (u, t) is set to be for m > Mj + 1. As we shall see in section 5, the r-functions 
of our principal interest (master T-operators for spin chains with rational i?-matrices) 
are, by construction, polynomials in u. This means that the coefficients aj m (u, t) are 
rational functions of u, so we assume this in the present section. 
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4.1.2 Krichever's conditions from the bilinear identity 



Since the poles at Pi are the only singularities of the product i/j u (t, z)if)*,(t', z) outside 
the contour C[o )00 ], the bilinear identity (3.24) then implies that 



N 



^res 2=Pi (V>„(t, z)if u ,(t' , z)) = 



i=l 



for all u,tj,u',t'j. Let us rewrite this identically in the form 



N 



J2pI u ' Ri(u',t';u,t) = 0, 



(4.4) 



i=i 



where the functions Rj are 



Mi 



m=0 



z=Pi 



From the polynomiality of the r-function it is clear that they are rational functions of vl . 



If so, the only way to satisfy (4.4) for all vl is to put each term equal to zero separately 



(because it is a sum of different exponential functions multiplied by rational functions 
which can be identically zero if and only if all rational functions vanish identically). 
Therefore, the bilinear identity is equivalent to iV conditions 



res z=Pi (ip u (t, z)^*/(t', z)) = 0, z = 1, 



or, more explicitly, 



Mi 

Y,~a<i m {u\t')d™{z- U 'e-W^ u {t, 



m=0 



1,...,N, 



(4.5) 



(4.6) 



Z=Pi 



which hold for any values of u, u' and t^t\. In particular, at vl = t\ = we have 

Mi 

= or 

z=Pi 



J]a m (0,0)9 2 m ^(t, 



z 



m=0 



Mi 



^2 a im d™ip u (t,z) 



m=Q 



(4.7) 



z=Pi 



where aj m = fcjai m (0,0) with arbitrary non-zero ki are parameters of the solution in the 



Krichever's approach. We call equations (4.7) the Krichever's conditions |17j . 
Remark. In the Krichever's approach, the parameters cii m can be taken arbitrary while 



in our case Oj m (0, 0) appear to be constrained by iV conditions of the form (4.23), as we 
shall see below. That is why we prefer to work with the set of independent parameters aj m . 
The initial values of the adjoint BA function coefficients 5i m (0, 0) are not independent 
and differ from the Oj m by the properly chosen N factors fcj. 



The set of conditions (4.7) yields a system of iV linear equations for iV coefficients 



Wk which demonstrates the consistency of the bilinear identity for BA functions of the 
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assumed form and allows one to fix them. Then the r-function associated with the ip- 



function according to (3.22) solves the MKP hierarchy. The points and the entries 



of the matrix cij m are free parameters of the solution. The coefficients w k appear to be 
rational functions of their arguments u, U while the r-function is a polynomial (multiplied 
by the exponential function of a linear combination of times). From the algebro-geometric 
point of view this solution is associated with a highly singular algebraic curve which is 
the Riemann sphere with cusp singularities at the points p^. 



4.1.3 Solving the linear system 



It is easy to see that conditions (4.7) are equivalent to the system of linear equations 

Ai(u,t) + ^2Ai(u-k,t)w k = 0, (4i 

k=l 

where 



M, 



A i (u,t) = J2^mdT{^ it ' z} ) 



m=0 



(4.9) 



Z=Pi 



These functions are polynomials in u, ti multiplied by the exponential factor p^e^ t,Pi \ In 
the case when all t t are equal to 0, they admit a more explicit representation: 

Ai(u,0) = Pi(a i0 + aap^u + a i2 pj 2 u{u-l) + a i3 pr 3 u(u-l)(u-2) + ...). (4.10) 

Note that ^(0,0) = a i0 . 



The system (4.8) can be solved by applying the Cramer's rule. This results in the 



following determinant representation for the BA function 

1 z~ l 
Ai(w,t) Ax(u-l,t) . 



-N 



. Ax(u-N,t) 



n ,, t; . A N (u,t) A N (u-l,t) ... A N (u-N,t) 



A^u -l,t) 



Ai(u - N,t) 



A N (u -l,t) ... A N (u-N,t) 



(4.11) 



From the definition (4.9) we have 



M k , 

A k (u,t- [z- 1 ]) = ^a fcm 9 c m Ce C(t,0 (l 



m=0 



A k (u,t)-A k (u + l,t)z \ 



C=Pi 



(4.12) 

Using this, it is straightforward to verify that equation (4.11) agrees with the general 



relation (3.22), with the r-function being given by the determinant in the denominator: 

At(u -l,t) ... Ai(«-JV,t) 



T u (t) = det Ai(u-j, t) 

i,j=l,...,N 



A N (u -l,t) ... A N (u-N,t 



(4.13) 
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N N 

It is a polynomial in u of degree ^ Mj multiplied by Y\ p^e^ t ' Pi \ 

3=1 i=l 



It follows from (4.11) that the last coefficient in (4.2), wn, is given by 



w N (u,t) = (-1) 



jyTit+l(t) 

rJt) 



We also note the property 



d tl Ai{u,t) = Ai{u + l,t) 



from which one can see that the first coefficient in (4.2), w\, is given by 



Wi(u,t) = -d tl logr„(t). 



(4.14) 
(4.15) 
(4.16) 



4.1.4 The adjoint BA function 



Similarly to (4.12), we have: 



A k (u,t + [z- 1 ])= J2a km d?(± 



m=0 



-C/z 



(4.17) 



Note that this function regarded as a function of z has a pole of order M k + 1 at z = p k . 
The principal term is 

M k \a kMkP u k +1 e^ 



A k (u,t + [z x ]) 



{z - Pk) 



+ 0((z- Pk )- Mk ) 



(4.18) 



We also see from (4.17) that the function A k (u,t + [z : ]) has a simple zero at z — 0. 

In order to obtain a more detailed information about singul aritie s of this function, 
let us add and subtract the term z u e^ t ' z ' ) in the numerator in (4.17) and separate the 
nonsingular part as z — > p k from the singular one: 

Z ( ^ a kmd( 



A k (u,t + [z- 1 ]) 



z 



u+lMt,z) 



E IIL 
(7 - 



ml a km 



m=0 



(z-p k ) m + l 



m=0 



z-C 



C=Pk 
(4.19) 

The first sum gives the multiple pole structure at the point p k while the second term is 
obviously regular at p k and has possible (essential) singularities and branching only at 
and 00. Note that A k (u,t + [z^ 1 ]) is, by construction, a rational function of z with 



simple zero at z = 0. This is obvious from (4.17) but not from (4.19). On the other 



which is rather obscure in (4.17). 



hand, the representation (4.19) explicitly shows the multi-pole structure of this function 



Rewriting (4.12) in the form 



A k {u, t + [z- 1 ]) = A k (u, t) + z~ x A k {u + 1, t + [z' 1 ]), 
it is straightforward to check that 



T u (t+[ Z - 1 ]) 



A 1 (u-l,t+[z- 1 }) A 1 (u-2,t) ... Ai{u-N,t) 
A 2 {u-l,t+[z~ 1 )) A 2 {u-2,t) ... A 2 (u—N, t) 

A 7V (u-l,t+[z- 1 ]) A N {u-2,t) ... A N {u-N,t) 



(4.20) 
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Expanding (4.17) in powers of z, we get 



A k (u, t + [z~ 1 ]) = A k (u, t) + A k {u + 1, t)z~ l + A k {u + 2, t)z~ 2 + ... 



and so the expansion of r u (t + [z x ]) around oo reads 



Ai(u+s-l,t) Ai(u-2,t) 
A 2 (u+s-l,t) A 2 (u-2,t) 

Ajv(w+S-l,t) Ajv(«-2,i) 



r„(t + K 1 ]) 



oo 

£ 

s=0 



2 



Ai(u-JV,t) 
A 2 (u-JV,t) 



(4.21) 



We thus see that the adjoint BA function has the determinant representation 

A l (u-l,t+[z~ 1 }) Ai(«-2,t) ... Ai(«-iV,t) 
A 2 (u-l,t+[z- 1 }) A 2 (u-2,t) ... A 2 (u-A,t) 



VC(M) 



z e 



-S(t,2) 



A JV (w-l,t+[z- 1 ]) Aiv(w-2,t) ... Ajv(w-iV,t) 



Ai(u-l,t) Ai(u-2,t) 
A 2 («-l,t) A 2 (u-2,t) 

An(u-IA) A N (u-2,t) 



Ai(u-N,t) 
A 2 (u-N,t) 

A N (u-N,t) 



(4.22) 



It has multiple poles at the singular points Pi of the algebraic curve which is the Riemann 
sphere with complicated cusp-like singularities. The principal parts at the poles are easily 



extracted from the determinant representation and equation (4.19) 



Since the function A k (u,t + [z *]) has a simple zero at z — 0, it is clear from (4.13) 
that the function r u (t + [z~ 1 ]) and thus the function z u e^ t ' z 'il>^ l (t, z) have zero of order 



iV at z = 0, whence the coefficients ai m (u,t) in (4.3) appear to be constrained by iV 
relations 



N Mi 

$^(-l) m (™ + n)\ a im {u, t)pl 

i=l m=0 



m— n— 1 



O n ,0 



71 = 0,1, 



N- 1. 



(4.23) 



4.1.5 The multi-pole structure of the adjoint BA function 

Let us introduce the notation 



A k (u,t) := det Ai(u + l-j,t) 

i — l, ...,]rf,...,N 

j=l,...,N-l 



(4.24) 



for the minor M k ^ of the NxN matrix Ai(u + 1 — j), 1 < i,j < N. Then, expanding 
the determinant in the numerator of (4.22) in the first column, we obtain: 

2-u e -e(M) N 



rjt) 



^(-l)*- 1 ^^ - 2, t) A k (u - 1, t + [z- 1 ]). 



k=l 
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Using (4.19) we can extract the poles: 



:(t, z) = ^(-i)*" 1 (M J' J E r _ °wi + terms regular at a11 (4 - 25) 

k=l m=0 



Therefore, 

^z= Pk [{z-p k ) n ^* u {t,z)\ = 
for all n > 0, or, in terms of the r-function, 



' rJt) 



(4.26) 



res 



[(z - Pk ) n z- u e-^T u (t + [z" 1 ])] = (-l)*" 1 ^ a kn A k (u - 2, t). (4.27) 



Note that the expression in the r.h.s. has the same structure as (4.13) and, therefore, is 
a polynomial r-function. This fact will be used below for the Backlund transformations. 

We also note the factorization formula 



res 



z=Vk 



(^(t'^)) = (-l)*- 1 ^^ A k (u, t) 



(4.28) 



which follows from (4.25). At v! — u — t'j — tj = this formula gives the relation 
between a k Q = afeo(0, 0) and a k Q. 



, ^ ^(-2,0) 
. — t) — — OfcO- 



ro(0) 



(4.29) 



4.2 Undressing Backlund transformations for the rational solu- 
tions 

As it will be clear in section 5, the main relations of the Bethe ansatz method are naturally 
built into the construction of rational solutions to the MKP hierarchy. In particular, the 
functions A k (u,t) and A k (u,t) will be identified, up to some irrelevant factors, with the 
(eigenvalues of) T-operators on the first and the last levels of nesting in the nested Bethe 
ansatz scheme. The parameters Pi will be identified with eigenvalues of the twist matrix 
g. We will also see that all these formulas can be understood in the operatorial sense (in 
the quantum space of the spin chain) since they involve only commuting T-operators. 

The nested Bethe ansatz scheme itself appears to be equivalent to a chain of some 
special Backlund transformations of the initial rational MKP solution with Krichever's 
data pi, a im that "undress" it to the trivial solution by reducing the number of singular 
points in succession. Here we present the idea solely in terms of the MKP hierarchy 
leaving the precise identification with objects from quantum integrable spin chains for 
the next section. 

Adding/removing a point pi to/from the Krichever's data of a rational solution is 
a Backlund transformation. It sends a polynomial r-function to another one. The key 



equation that allows one to implement such a transformation is (4.27). Let us start with 
the n = case. Specifically, consider the function 

rM(t) = (-lrWp, (z- u ~ x e-^ z h u+1 {t + [z- 1 ])) . (4.30) 
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According to the n = case of (4.27), it is equal to 



r!Ht) 



Ai(u-l,t) A 1 (u-2,t) 



Ai-x (u — l,t) 



.4 



t+i 



Ut-l,t 



A N (u-l,t) 



A i _ 1 («-2,t) 
A i+1 («-2,t) 

Ajv(tt-2,t) 



Ai(u-JV+l,t) 

A_!(u-iV+l,t) 

A m (w-iV+l,t) 
A N (u-JV+l,t) 



(4.31) 



i.e., to the minor of the matrix Ai(u + 1 — j) multiplied by ajo- Therefore, it is a 
r- function, i.e., it satisfies the same 3-term Hirota equation as r u (t) does and r — > is 
indeed a Backlund transformation. Below we assume that a i0 7^ for alii = 1, . . . , iV, 
then the transformation is non-trivial. 



Having at hand t« (t), one can construct the BA function 

^(M) = ^ ^ 1(t - [ f D 

Tiz ( t ) 



(4.32) 



which obeys the same Krichever's conditions (4.7) at the points pi,. . . ,p^ except the 



point pi, where no condition is imposed. Note that the BA function ijM has the expansion 



around the point z = of the same form (4.2) but with a pole of order N — 1 rather 



than N, so the number of conditions again matches the number of unknown coefficients 
of the singular part at z — 0. The adjoint BA function is 



^T(t } z) = z- U e 



«„-€(M) ^(t + t*- 1 ]) 



r£ J (t) 



(4.33) 



It has multiple poles at the points pi, . . . ,Pn except the point pi. 

This process can be continued by taking the residue of z _w-1 e - ^ t,2 V^ 1 (t + [z^ 1 ]) at 
some other singular point pj, j =fi i: 



-l)^res 2=p . [z- u ~ l e-^T [ : ] +1 (t + [z- 1 ]) 



(-ly^-^res^ (z^y^e-^-^izj-zAr^it + [ Z ; 1 ] + [zj 1 ]) 



z j=PJ 



(4.34) 



aiofljo det A r (u — s,t) 

l<r<N,^i,j L 
s=l,...,N-2 



Here and below e^- = 1 if i < j and —1 if i > j. 

We thus obtain a chain of Backlund transformations r — > r' 11 ' — > r^ 1 * 2 ' — > 
ereneral recursive formula at the n-th level is 



.. The 



[ii...i r , 



\ in — 1— Card{fc|ife<i n ,l</c<rt— 1} 



res. 



(z~ u -\ 



e(t,z in ) T [»l-jn-l] 



u+1 



(4.35) 
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which leads to 

T [h-in]( t ) = (_l)^ rrs 



Z 'k~ Pi k 
Kk<n 



n 

(n 

a=l 



n 



a i a 



n 

jA n (z h ,. . . , Z in )T u+n (t + J^k^]) 



) det 


A r (u — s, t) 


1 l<r<N,^h,...,in 




s=l,...,N-n 





In the first line d n = i± + . . . + i n + |^(^ + 1) and 

A n (z h ,...,z in ) = \ [ 



(4.36) 
(4.37) 



i>j 
i,j£{ii,~.,i n } 



is the Vandermonde determinant. In particular, on the second to last level we have 



N 

t$* -^(t) = ( J] a^Mu -l,t) 



(4.38) 



and on the last level rl 12 "'^(t) = flili a «o = const. Note that a i0 = Ai(0, 0). 



By taking residues of the Hirota equations in the form (3.7), (3.9), one can obtain 
bilinear relations for the r-functions on the neighboring levels of this chain. In particular, 
we have: 

^(^(t) + e jk r^\t)rf +1 {t) + e K fl(t)4(t) = 0, (4.39) 
s^(t)r u+1 (t) = ^(t^Ut) - T^(t)r£Ut). (4.40) 



We can also rewrite the determinant (4.36) in terms of (4.38): 

2V 



T, 



il...i r , 



>= (n 



-N+n+l 



det 

l<r<N,^ii,...,i n 
s=l,...,N-n 



(4.41) 



Note that (4.39) and (4.40) are respectively the Pliicker and Jacobi identities for (4.41). 

Equation (4.27) allows one to make the same Backhand transformations by picking the 
coefficients in front of higher order poles of the function r u (t + [z^ 1 ]) at z = pi, . . . ,p^- 
The results differ by normalization factors independent of u, t. For example, instead of 
taking residues one may pick up the highest singularity (the pole of order Mj + 1 at Pi) 
and define 



r«(t) = (-l^-ires^f^-pi)^ 



\Mi 7 -u-l p -£{t,z) 



r u+ i(t + [z- 1 ]) 



(_l)«-ip-«- 1 e-«*»w) res z=Pi [{z- Pi ) M H u+x {t + [z' 1 }) 
a iMi Ai{u - l,t). 



(4.42) 
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In a similar way, one may define the "undressing chain" of such transformations as in 
dOB] and fl436l): 



T (ii...i n ) _ 1— Caxd{fe|ife<i»,l<fc<ri— 1} 



x res,. — 7v 



- P J M ^zC~ 1 e-<^T^ l) (t + K 1 ])) , (4-43) 
which leads to 

n 

rt- in \t) = (J] MJ a iaMiQ ) det \A r (u - s, t)l . (4.44) 

s=l,...,JV-n 

The construction of the undressing based on the highest poles is somewhat more compli- 
cated because it explicitly uses the multiplicities of the highest poles but the advantage 
is that it always gives the non-vanishing result (because the coefficients in front of the 
highest poles cannot vanish by their definition). 

4.3 Fermionic realization of rational solutions 

Here we suggest a realization of the rational solutions as vacuum expectation values of 
some fermionic operators. 

As before, we fix N points pi and the coefficients Oj m . Let us introduce the following 
fermionic operators: 

Mi 

*i{z) :=^a im cT#z) (4-45) 

m=0 



The operators ^i(pi) are particular cases of (3.11) when each set contains just one 



element N + 1 — i. Therefore, we already know that the matrix element 

r„(t) = (n\ e J +^ 1 (p 1 ) . . . * N (p N )\n-N) (4.46) 

is a r-function. 



On the other hand, we can calculate it using the Wick's theorem in the form (2.26). 
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We get: 



T n (t) = det (n-j + l\e J +®%(pi)\n-j) 

i,j=l,...,N 



■ det M J2 a ™ d ™ \- e({t ' Z) (n-j+l\^(zWn-j \n-j+l)] 

i,j=l,... ,N *— • 
m=0 



m=0 fcSZ 



(4.47) 



'J 



m=0 



det Ai(n-j,t) 

i,j=l,...,N 



which coincides with the r-function (4.13) after the change n — >■ u. Let us also note the 
formulas 



T, 



n (t-lz' 1 ]) = z - n e~^ (n + l\ e J +«V(^i(Pi) . . . * N (p N )\n-N) 

r n (t+iz- 1 ]) = z n ~ l e^ (n - 1| e J +<V(z)*i(Pi) • • • **(p*)| n-iV) 
which directly follow from the bosonization rules (3.4). 

The simplest example is Mj = when ^i(pi) = tpipi) (we take a i0 = 1) and 

r (M l= o) (t) = ( ri | e ^(t)^( Pl )...^( pjv )| n _iV) 

N N 



(4.48) 
(4.49) 



which is easily calculated to be r { n M ^\t) = JJ^" 1 ~ Pj 1 ) \\pl~ 1 e i{t ' Pk 



i>j 



k=l 



In the same way as (4.46), we have a fermionic realization of the intermediate r- 
functions for the Backlund transformations. The r-functions corresponding to (4.36) 
have the form 



rt- lm] (t) = (J[a ia0 ){n\e j +V JJ y^n-N + m) 



(4.50) 



a=l 



j=l,^ii,...,ir, 



where m is a non-negative integer (0 < m < N). The r-functions corresponding to (4.43) 
have the same form except for the prefactor (■■■). These r-functions reduce to (4.46) 
for m = 0. The determinant formulas (4.36) for these r-functions follow from the Wick's 
theorem (2.26) in the same way as m = case (4.47). The BA functions for (4.50) are 
defined by 



(t,z) = z n e 



nMt,z) T _n 



r, 



- Am] (t) 



(4.51) 
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The BA functions for (4.43) are defined exactly in the same way. In addition, both of 
them coincide under the condition that the prefactor of the r-functions are not zero. By 
using a formula similar to (4.48), we obtain 



(n + 1| e J +<*ty(*) Uf=i,^ u ..., im *i(P/)l n-N + m) 



l\...l r , 



(n| e^(*) njli,^,...,^ *i(Pi)| n-JV + m) 
Then the Krichever's conditions for the Backlund transformations 

for £;G {l,2,...,N}\{h,...,i m } 



j=0 



z=Pk 



(4.52) 



(4.53) 



follow straightforwardly from (4.52) and the relation (^k(Pk)) = 0- 



5 The r-function in quantum integrable models 

5.1 T-operators 

We consider generalized quantum integrable spin chains with rational GL(iV)-invariant 
i?-matrix. The i?-matrix acting in the tensor product of the iV-dimensional vector rep- 
resentation of GL(N) and an arbitrary finite dimensional representation 7r\ (labeled by 
a Young diagram A) has the form 

N 

R x (u) =«I<g>I+ eji®7r A (ey), (5.1) 



where u G C is the spectral parameteijj and are standard generators of the algebra 
gl(N) with matrix elements (ejj) a/ g = 5i a bjp- In the second term, the vector representa- 
tion (the one for the Young diagram with one box) 7rn(ejj) is abbreviated to eji. 

A large family of commuting operators (quantum transfer matrices or simply T- 
operators) can be constructed as 

T\u) = tr^ (R x (u - ® • • • ® R\u - ^)(r L ® 7r x (g))) , (5.2) 

where ^ are arbitrary parameters (inhomogeneities at the lattice sites) and g G GL(N) is 
called the twist matrix. For technical simplicity, we assume below that the twist matrix 
has distinct, non-zero eigenvalues. The trace is taken in the auxiliary space where the 
representation tt\ acts. The T-operators act in the physical Hilbert space of the model 
= (C N )® L . The Yang-Baxter equation implies that the T-operators with the same g 
commute for all u and A and can be diagonalized simultaneously. It follows then from 



(5.1) and (5.2) that all their eigenvalues are polynomials in u of degree < L. 

It is natural to define the T-operator for the trivial representation (empty Young 
diagram) as 

L 

T\u) = H(u-^), (5.3) 

i=l 



i Om definition of the spectral parameter differs from that in H] by a factor of 2 
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where the multiplication by the identity operator in the r.h.s. is implicit. We introduce 
a special notation 

L 

0(u):= Y[(u-£ t ) (5.4) 

i=l 

for the polynomial in the r.h.s. This polynomial characterizes the model through its zeros 
£j. Another normalization, sometimes more convenient, is to set 

TV) = £g. (5.5, 

In this normalization all T-operators are rational functions which may have poles only 
at the points & and T (w) = 1. 



5.2 Determinant formulas for T-operators 

The T-operators for different representations in the auxiliary space are known to obey 
an infinite number of functional relations which follow from the fusion procedure and, 
eventually, from the Yang-Baxter equation. 

First, there are determinant formulas which express T x (u) for arbitrary A through 
T-operators for symmetric representations T s (u) := T^(u) corresponding to one-row 
diagrams (s) [HUE]- They are called Cherednik-Bazhanov-Reshetikhin (CBR) formulas 
or quantum Jacobi-Trudi identities. In our normalization they have the form 

/AH \ ^ 




T\u) = [ | 4>(u — k) det T Xi „ i+j (u + (5.6) 

ij=i,...,\ 1 

A direct proof "from the first principles" for the models with the rational GL(N\M)- 
invariant i?-matrix was given in [5]. The proof uses the operation of co-derivative with 
respect to the twist matrix g (see below). When the dependence on the spectral parameter 
disappears (say, in the limit u — > oo) one obtains the standard Jacobi-Trudi formulas 



(2.13) for characters. We note in passing that the determinant formula for the generalized 



Schur functions (constructed in terms of the Gauss decomposition of a matrix (resp. 



orthogonal polynomials)) obtained in [30] (resp. [31]) has a structure similar to (5.6). 

Second, there are determinant formulas of the Giambelli type which express T x (u) 
for arbitrary A through the T-operators T^('u) := T" +1,1 \u) corresponding to the hook 
diagrams (I + 1, l fc ) 



T\u) = (<f>(u)r diX)+1 . . det T Ai _. A ,_ i(u ). (5.7) 

i,j=l,...,d(\) 3 

Note that there are no shifts of the spectral parameter here. In Appendix C we prove 



that the determinant representation of the Giambelli type (5.7) follows from the quantum 



Jacobi-Trudi identities (5.6). 



In the normalization (5.5) the pre-factors cancel and the determinant formulas look 
simpler: 

J x (u) = det T\._ i+ j{u — j + 1) = det J x ^ y (u). (5.8) 

i,j=l,...,X[ i,3=l,...,d(X) 3 
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Below we will argue that these formulas allow one to establish a close connection with 
classical MKP integrable hierarchy and to recover a hidden free fermionic structure in 
the auxiliary space. 



5.3 The master T-operator 



Let us introduce a generating function of the T-operators (the master T-operator) de- 
pending on an infinite number of parameters t = {ti,t%, . . 



T(u,t) = J2sx(t)T> 



u 



(5.9) 



These operators commute for different values of the parameters: [T(u,t), T(u', t')] = 0. 
The T-operators T x (u) can be restored from the master T-operator according to the 
formula 

T\u) = s x (d)T(u,t) 



(5.10) 



t=o 



where d = {d tl , \d t2) \d t% , . . . }. In particular, T(u, + [z x ]) is the generating series for 
T-operators corresponding to the symmetric representations: 



T(u, [z- 1 ]) = 5>-*T s 



u . 



(5.11) 



s>0 



5.3.1 The master T-operator as a r-function 



Comparing formulas (5.6) with (3.18) and (5.7) with (3.17), one can see that they coincide 



after the identification u «-» n and T (u) -H- C\(n). According to section 3.3.2, this implies 



that there exists a fermionic group-like element G (an operator in the quantum space of 
the model) and a ^-dependent operator C{u) such that 



T\u) = (-l) bW C(u)(X,u\G\u) 



(5.12) 



(which makes sense at least for integer values of u), and thus the master T-operator is 
the r-function of the MKP hierarchy: 



T(u,t) = C(u) (u\e J +^G\u) 
Setting t k = 0, we get T(u, 0) = T (u) whence 

C(u) (u\G\u) = <j){u) 



(5.13) 



(5.14) 



(multiplication by the identity operator in the quantum space of the model is implied in 
the r.h.s.) 

A weaker statement, that the master T-operator is the r-function of the KP-hierarchy 
if T x (u) in (5.9) has the CBR-determinant representation, can be proved independently 



in the following way. Suppose T (u) has the CBR determinant representation (5.6). 



Then it can be expressed as the Giambelli determinant (5.7) in the Frobenius notation 
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for the Young diagram A = . . . , a^|/Si, • • • , A*) (see Appendix C). For any k, I G 
{1, 2, . . . , d(\)}, k < I, the Jacobi identity for this Giambelli determinant reads: 



rp(ai,...,& k ,...,&i,...,a d \Pi,... ,/8fe,... ,A,...,/3 ti )^ M -j2 n (ai,...,a d |/3i,...,/3d)^ M ^ 

_ J>(ai,...,afc,...,ad[j8ir" A,-,fid) (^A'J'icti ,—,&l,—,<Xd\0i,—,Ph —M (y\ 

_|_ 2i(«i,...,a fc ,...,a d |^i,... ->Pd) ^ u ^Ji(ai, —Ah— ,<Xd\fh>-»tPkr"tM fy\ = Q (5 ]_5^ 

where (resp. means that at (resp. /?/) is removed from the original diagram 
(ai, . . . , CKd I /?i , . . . , (3d). This is nothing but the 3-term relation given in [II] . which is 
known to be the necessary and sufficient condition that T(u, t) solves the KP hierarch 
Therefore, T(u, t) is a r-function. In other words, the functional relations ( 5.6[ ), ( |5.7| 



5 



mean that the T-operators T x (u) are Plucker coordinates for the r-function T(u,t) 
We emphasize that the notion of the master T-operator is independent of the specific 



form of the i?-matrix[j The key point is that T x (u) in (5.9 ) is the transfer matrix obtained 
by the standard fusion procedure in the auxiliary space (the horizontal leg) labeled by 
the Young diagram A. The master T-operator is thus well-defined for spin chains with 
trigonometric and elliptic i?-matrices as well and the statement that it is a r-function 
still holds. The statement is also applicable to the T-operator corresponding to any 
higher representation in the quantum space of the model. It can be also applicable to 
the g-characters of representations of quantum affine algebras [35] . 



5.3.2 Hirota equations for the master T-operator 

As any r-function, the master T-operator satisfies the bilinear identity 

j> e^'^z^'Tiu, t - [z- 1 ]) T(u', t' + [z~ 1 ]) dz = 



(5.16) 



for all u, u',t,t'. Here the integrand should be interpreted as a formal power series on z 
and the contour C should encircle z = oo. By standard manipulations, one can derive 



from (5.16) the infinite KP and MKP hierarchies of differential equations. Making the 
shift of variables t — > t + a, t' — > t — a in (5.16), where a = (a 1; a2, . . . ) and setting 



5 at and /3fc here are related to the parameters in page 959 of [16. as ik = otk>jk — —fik — l,r = d. 
6 For the trigonometric case (or the rational case) , we can define the "universal master T-operator" . 
Let R be the universal R-matrix associated with U q (gl(N)). R is considered to be an element of B + ®B_, 
where B± are Borel subalgebras. Let tt\(u) be an evaluation representation of U q (gl(N)) based on a 
tensor representation of U q (gl{N)) labeled by the Young diagram A. This representation depends on 
the spectral parameter u since the evaluation map does. We can define the universal T-operator (cf. 
[33] ) as T x (u) = (Tr ® id)(n\(u) <g> id)((g <8> 1)R), where g is a twist element (a function of the Cartan 
elements). Then our universal master T-operator is obtained by substituting this into (5.9). Now (5.9) is 
an element of £?_ . One can obtain a master T-operator on a particular quantum space V if we consider 
a representation of B-. For example, if we consider V = 7T(i)(£i) ® tt(i)(£2) ® ••• ® tT(i)(£l), where 
{tt(i) (^ j)}^— i are fundamental representations we obtain the master T-operator for a lattice model in 



section [5.3.3 (as a limit q 1). If we consider a space of the vertex operators of CFT for V, we will 
obtain a master T-operator for CFT. It is tempting to rewrite the universal master T-operator as a kind 
of determinant over some function of (tt^(u) ® id)(i2), where nru(u) is the fundamental representation. 
This should be a generalization of a generating function of transfer matrices (cf. pM]L 
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u — v! G Z> , one obtains, in the usual way, the following Hirota-type bilinear differential 
equation for the MKP hierarchy: 



h m (2eL)h m+u _ u , +l {-i>)e^ akDt KT{u,t) ■ T(u',t) = 0. 



(5.17) 



m=0 



Here TJ t = (D tl , \D t2) \D t . A) . . .), where D tk is the Hirota derivative with respect to t k . 
For example, the first equation of the MKP hierarchy reads 

T(u + l,t) , . , , x \ . (n , T(u + l,t) 



T(u,t) 



The bilinear relations (3.7), (3.9) for T(u,t) can be written as 

(z 2 - z 3 )T (u,t + fe- 1 ]) T («,t + fa 1 ] + [23- 1 ]) + (231) + (312) = 



(5.18) 



(5.19) 



z 2 T (u + 1, t + [zf 1 ]) T (it, t + [z 2 1 ]) - z x T (u + 1, t + [zj 1 ]) T (m, t + [zf 1 ]) 

+ (z 1 -2 2 )r(M + i,t + [^ 1 - 1 ] + [2 2 - 1 ])r^,t) =0 



(5.20) 



(comparing to (3.7), (3.9), we have made the overall shift of the time variables t — > 
* + \ z i l \ + feT 1 ] + [«s X ] in d 5 19 l ) and t ^ t + f^ 1 ] + [z^ 1 ] in fl5 20| )). 

For any positive integer if and any subset i 2 , ■ ■ ■ , C {1,2,..., K} where all ik 
are different, we introduce the operator 



T{M 2 ,..M( M )=TM + XK 1 ])- 



fc=i 



Solving (5.20) for (5.21) recursively, we obtain the determinant formulas 

deti< kJ < n (zj; n T^{u-j + l,t)) 



rp{l,2 > ...,K}\{%un,...,in} ■(-) 



deti< fcjj < n {z\ h n ) ]1Li T ( m - M) ' 

detx< fcj -< w (z{; n T^-> K ^(u + j - l,t)) 
det x <^ B {zf; n ) UVi Vi (« + fc, t) ' 



(5.21) 



(5.22) 



(5.23) 



(5.22) 



Note that (5.19) is a Plucker identity for (5.22) and (5.20) is the Jacobi identity for 



5.3.3 The master T-operator in models with rational i?-matrices 



For models with rational -R-matrices, the T-operators T (u) (5.2) admit a more explicit 



realization in terms of the co-derivative operators acting on functions of the twist matrix 
g G GL(N) [9]: 



T\u) = (u-fr + D) 



{u-i L + b) X x{g) 



{u-& + D) X x{g), 



(5.24) 



8=1 
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where X\{d) — ^ I7 ^x{g) is the character of g in the representation 7r A . All necessary 
information about the co- derivative D is given in Appendix D. 

Set 



F fl (t) = exp(^* fc tr</*) . 

\fc>i / 



(5.25) 



As immediately follows from its definition, the master T-operator (5.9) can be written 
in terms of the co-derivative as follows 

L 

T(u, t) = (g) (u - & + D)F g (t). (5.26) 

i=l 

Note that each term in the direct product introduces a spin into the underlying quantum 
spin chain. 

Alternatively, using the fact that (det g)~ u D(det g) u = D + u, we can also write it as 

L 

T(u, t) = (det g)- u (g) (D - &F g (u, t), (5.27) 

i=l 

hiding the spectral parameter into the trivial r-function (with no spins) 
F g (u,t) = (detg) u F g (t) = (det gf exp t k tr g k ) . 

k>l 

In this form it is clear that the spectral parameter is naturally combined with £j into a 
hierarchy of times u — to, ti, i2> • • - : 



F g (u,t)=exp(j2tk^9 (k) ), where = { 



g fe if /c > 
log^ if = 0. 



(5.26 



fc>0 

The explicit form of the simplest T-operators for a chain of length L = 1, 2 is: 



T (L=1) (u,t) = ( (u 



- eoi + X) J ex P (X t; tr 9 1 ) > ( 5 - 29 ) 

fe>l / Z>1 

T(l =2) (m, t) = ( (u - e0(w - 6) I ® I + (u - £i) X I®g k + (u-&)J2 kt k g k ® I 

V fc>i fc>i 

fc-1 \ 

+ E E ^ ® ^ + E E Hfc p ^ ® s*~ a ) ex p (E *' tr ^) > 

fc>l fc'>l fe>l a=0 / Z>1 

(5.30) 

where P is the permutation operator (see Appendix D for definitions). For clarity, the 
unit matrix I is everywhere written explicitly. 

Using the method of [8], it is possible to show that the master T-operator contains 
the whole family of Baxter Q-operators as well as the T-operators that emerge on inter- 



mediate levels of the nested Bethe ansatz. This will be explained in section 5.5 
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It follows from (5.26) that 



T(u,t-{z- 1 }) = [£){v 

i=\ 



-Zi+D) [(det(I-2-^))F 9 (t 



(5.31) 



T(u,t+ [z- 1 ]) = (g) (u-&+ D) [(det(I - z' 1 g))' 1 F g {t)] . 
i=i 

Using these formulas, one can explicitly verify that 

lim U-z) ±N T(u } t t [z- 1 ])} = (det g^Tfa ± 1, t). 

z— >0 



(5.32) 



Comparing with the prescription (3.14), we see that there is an extra factor (detg) 



which is due to a slightly different definition of the r-function: in our case it is a pure 



polynomial in u while in (3.14) it is a polynomial multiplied by an exponential function. 



Denoting f) L {u) := (£)f =1 (u — D) and w(z) = (det (I — zg)) 1 for brevity, we can 
rewrite (5.19), (5.20) in the form 



(z 2 - z 3 ) (v L (u)w(z^)F^ (v L (u)w(z^)w(z^)F^ + (231) + (312) = 0, (5.33) 



z 2 [V L (u + l)w(z^)F g ) (VLiuMz^Fgj-zAVLiu + lMz^FMVMwiz^F, 



+ { Zl - z 2 ) ( V L [u + l)^^- 1 )^^- 1 )^ ) ( V L (u)F g 



0. 



(5-34) 

In the second equation, one can recognize the "master identity" from [8] (to identify 



them, the redefinition zi — >■ Iffy is required). The relation (5.32) implies that equation 



(5.34), which is the Hirota equation for MKP hierarchy, follows from (5.33), which is the 



£l = 0, the third term in equation (5.34) vanishes since 



Hirota equation for the KP hierarchy, in the limit z^ — > 0. Note that at t{ = (i.e., 
F g = 1) and u = £i = • • 

T>l -1 = and the equation becomes the "basic identity" (5.9) from [U] 

z ((l + D\® L w(z) \ [b® L w{z') \ = z' ((l + D\® L w(z')) (f)® L w(z)) . (5.35) 



where we put z 



z 1 ,z 



When z — > 0, (5.35) reduces to the following identity [9] linear in w, 



D® L w(z') = z' {gxV^L) Ul + by L w{z') 



where we used w(z) = 1 + zti g + 0(z 2 ) and 

t)® L tr g = g{P x ... L . 



(5.36) 



(5.37) 



Here V\,.x = V\ 2 V2-i • • • Vl-i,l is the full cyclic permutation in the quantum space stem- 



ming from (D3) and from the repeated use of |D7J) (see Appendix D). The subscript in 



gi means that it acts, as a matrix, in the space number 1 of the spin chain. Note that 
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the order of operators in (5.35) is irrelevant due to the commutation of transfer-matrices 



with arbitrary auxiliary representations. Similarly, for z' — > 0, in the first order in z' we 
obtain: 



z 1 



D 



w[z) 



D 6 



w[z, 



{gxPx...i 



(5.38) 



It is a direct consequence of (5.36) if we recall that the two factors in (5.38) commute: 



both are particular cases of the same family of commuting transfer-matrices. Eqs.(5.36) 



and (5.38) were proven in a combinatorial way in [9]. Multiplying the left hand sides 



and the right hand sides of (5.36) and (5.38) we restore again the relation (5.35) which 



served there as the main tool for the direct proof of the CBR formulas. 



us now show how the case F g = 1 of (5.34), i.e. 



As is explained above, equation (5.38) was proven in [5] in a combinatorial way. Let 



z 2 [f> L {u + lMzr 1 )) [V L {u)w{z 2 1 ^ - Zl (v l (u + l)w(z 2 l j) (t 



L (u)w(z l 



+ ( Zl - z 2 ) (f) L (u + lMzrX^ 1 )) {Vl{u)1) = 0. (5.39) 



can be deduced from it. First we notice that when L > 1, equation (5.35) is the particular 
case of (5.39) when £j = u for all % = 1, . . . L. Next we notice that although (5.35) only 



holds if L > 1, the third term in (5.39) is such that (5.39) also holds (trivially) when 
L = 0. Since the left hand side of (5.39) is polynomial in the variables Ui := u — £j, this 



allows for a recurrence procedure which proves that the term with an arbitrary degree in 
each Ui — (u — £j) in (5.39) is equal to zero. See Appendix E for a detailed proof on this. 



5.3.4 Proof of the determinant representation and the Hirota equation for 
the master T-operator 

A chain of arguments in (9j E] allows one to give a self-contained proof "from the first 
principles" of the central claim of this paper: the master T-operator for rational spin 
chains is a tau-function of the MKP hierarchy. Here we give a direct proof of this 
statement using the representation of the master T-operator through the co-derivatives. 



The master identity in the form (5.20) or (5.34) directly follows from the determinant 



formula (5.22) by applying the Jacobi identity. So in order to prove the master identity it 



is enough to prove (5.22). In fact, for this proof we can assume, without loss of generality, 
that t = 0. Indeed, the product^] 



K 
k=l 



K 



X\w{z k ) 



(5.40) 



t=o 



fc=i 



can approximate F g (t) for any set of t with any needed accuracy, if K is sufficiently large, 
by an appropriate choice of the parameters Z\, z 2 , ■ ■ ■ , zk- In other words, we change the 
variables from t to z\, z 2 , ■ ■ ■ , zk (which is the Miwa change of variables): 



k \ K 
T^ 2 - k \u):=t[u,Y,[^} \ = V l (u) ]Jw(z k ) 

k=l J k=l 



(5.41) 



7 Here we change all —t ~ for convenience. 



38 



or, for a subset {ii, i 2 , ■ ■ ■ , i n } C {1,2,..., K}, 



k=l 



V L (u) l[w(z ik ). 

k=l 



(5.42) 



Then, in order to prove (5.34) for arbitrary t, it is sufficient to prove (5.34) for arbitrary 
t = + Xl/b-i ' which is equivalent to proving (5.22) when t = 0. Therefore, we only 
have to prove the identity 



J-i{il,i2, —,in} ! 



u 



deti< fe)i <„ [Zi k 3 V L {u - j + 1) w{z. 



deti< fcJ < n {z\ 3 ) nLi Y\i=ii u -Zi-k) 



(5.43) 



This formula expresses the general T^ %1,l2 '"' ,%n ^{u) through TW(u) with only one index 
(the generating operator for transfer matrices in symmetric representations). The proof 
goes along the same lines as the proof of the CBR formula in |9J. It consists of two steps: 

1. We should prove that Tf^.-.M 

(u) is a polynomial operator of the power L, as it 

is clear from the definition; 

2. We should prove that the coefficients of polynomials in the r.h.s. and the l.h.s. of 
(5.43) are equal. 

In order to prove the first statement, it is enough to show that the determinant in 
the numerator vanishes at all zeros of the denominator, namely at u — fa + j, I = 



, n— 1. For that, let us compare two consecutive columns, the j'th 



and the (j + l)'th, of the matrix in the numerator. If the 2x2 minor determinant 



det 



z n ik 3 V L [u - j + 1) w(z ik ) z™ 3 V L (u- j) w{z ik ) 
zZ j T> L {u-j + l)w{z ikl ) zZr X V L 



[u- 3)w{z Kli 



(5.44) 



built out of any two arbitrary rows k and k' of these two columns are zero, then the 
whole determinant vanishes^} Indeed, the corresponding operatorial identity 



z (V L (b + l)w(z)) (V L (^)w(z')j - z' (V L (& + l)w(z')) (V L (^)w(z)j = (5.45) 

appeared in [U] in relation to the proof of the CBR formula. This identity can be deduced 
from (5.39), where the third term is equal to zero at position u = fa. The identity (5.39) 
itself is, as shown above, deduced from the relation ( 5.35[ ) (which is proven diagrammat- 
ically in [9]) by recurrence over the size L of the spin chain. 

The proof of the second statement, about the equality of coefficients in polynomials 
on both sides of (5.43), can be achieved by comparing their large Ui = u — fa asymptotics. 
First, we can rewrite (5.43) as follows: 



//j u-2 ' • • ul deti< fe j<„ ( zf 3 ®f =l ( 1 



i 



D w(zi 



det 



l<k ,j<n 



8 Of course these minors vanish only under the constraint that u = £i +j, where 1 = 1, 
j = 1, ...,n- 1. 



(5.46) 
, L, and 
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The r.h.s. must be a linear polynomial in each of the variables U[. We can then reconstruct 
this polynomial from the large ui asymptotics. For example, for large U\ we find 



u 2 u 3 ■ ■ ■ u L {ui + D) deti< fcii <„ izi k 



tf=2 



1 



1 



Ul+l-j' 



D)w(z 



det 



l<k,j<n 



(C) 



(5.47) 



where we used the fact that at large Mi the co-derivative in the bracket ( 1 H , , 

1 + ^Dj can occur only once in each of L\ terms of determinant and the Leibniz rule 

to give the right overall asymptotics U\ x l/u\ = 1. Expanding in this way for each of 
the remaining u^s we recover in the r.h.s. the original definition (5.42). 



In this way we have proved the determinant formula (5.22). It amounts not only to 
the proof of the master identity (5.20) but also to the fact that the master T-operator 
satisfies the whole general Hirota bilinear identity (5.16) for the MKP hierarchy. 



Let us note that the commutativity of two general master operators (1.5) follows, in 
virtue of the determinant representation (5.22), from the commutativity of generating 
functions of T-operators for the symmetric tensor representations 



V L {u)w{z)) [V L (u')w(z' 



V L (u')w(z')) [V L (u)w(z 



(5.48) 



A proof of this last relation can be also found in [9] . 



5.3.5 Analyticity properties of the BA functions related to the master T- 
operator 



Equations (5.31 ) also allow us to characterize the class of solutions of the classical Hirota 
equations relevant to the quantum spin chain. First of all, it is obvious that the master 
T-operator is a polynomial in u of degree L. Second, consider the (operator-valued) BA 
function and its adjoint: 



T(u,t) 



z e 



r u (t,z) 



z -u e -{(t,*) 



Tfat + jz- 1 }) 
T(u,t) 



(5.49) 



From (5.31) it is clear that the function ip is of the form (4.2) while the function ip* is 
of the form (4.3) with the poles at the points pi which are the eigenvalues of the twist 
matrix g (assumed to be all distinct and nonzero). This means that in the appropriate 
basis g can be diagonalized: 



g =diag(pi,p 2 , 



,Pn) 



(5.50) 



The maximal possible order of each pole is L+ 1. Indeed, the T-operator for the empty 
spin chain (at L = 0) gives first order poles at p^. 



N 

T (L=0) (u,t+[ Z - i ])= (n 

i=i 



Pi 



N 



i=l k>l 



(5.51) 



and each co-derivative increases the order of each pole by 1. Note also that (5.51) (and 
thus the adjoint BA function) has zero of order N at z = 0. 
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5.4 The wave operators as non-commutative generating series 

The wave operators (3.30), ( |3.31[ ) for the master T-operator are of the form 



W-\u,t) = ye-^ hs ^ U + \ t 



(5.52) 
(5.53) 



Expanding both sides of (5.31 ) in powers of z, we get 

L 

h k (d)T(u,t) = (9) (u-&+b) bck(g)F g (t)} , 



h k (-d)T(u,t) 



1=1 



II- 



[(-l) k x k (g)F g (t 



(5.54) 



(5.55) 



i=l 



where Xk{g), X (ff ) are . characters of symmetric and antisymmetric representations re- 
spectively. From (5.55) it is clear that h a (—d)T(u,t) = for a > N because the anti- 
symmetric GL(N) character x a (g) vanishes. Therefore, the series (5.52) is truncated: 

N 

1 1 i / M / / ii ± \ 

(5.56) 



v ' ; ^ T(u,t) 



a=0 



Moreover, since for GL(N) x (q) = det (7, we have from (5.55): 



h N {-d)T(u,t) 



[u- 



-Zi+D) [(-l^det^t 



i=i 



-l)"detgte)(u+l-Z i +D)F g (t) 



l) N det gT(u + l,t) 



and so the coefficient in front of the last term in (5.56 ) is equal to (— 1)^ det g - \ — 

1 I T(u,t) 

as it should, as soon as the points pi are identified with eigenvalues of the matrix g 
(compare with (4.14) and recall that the r-function from section 4 and the master T- 
operator from this section should be identified as r u (t) = (det g) u T(u, t)). 



We conclude that the r-function (5.9) corresponding to the GL(iV)-invariant spin 
chain on L sites must obey the following two conditions: 

• It is a polynomial in u of degree L multiplied by exponential function of a linear 
function of times, 

• The series for the wave operator (or, equivalently, the expansion of the Baker- 
Akhiezer function) truncates at the iV-th term. 
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At last, we remark that the wave operator (5.52) at £j = is nothing else than the 
"non-commutative generating series" of T-operators for antisymmetric (fundamental) 
representations T a {u) := T^ ia \u): 



N 



T a (u) 



w(u, o) : = w(u) = J2(-i) a -rr e ~ adu - 

<p(u) 



a=0 



(5.57) 



Similarly, W 1 (u) is the "non-commutative generating series"]^] of T-operators for sym- 
metric representations T s {u) := T^ s \u) (cf. [3]): 



w- i (u) = jre- sduTs(v 1 



s=0 



1 z - 



T s (u-s + l) 
(j){u-s+l) 



(5.56 



The operator T N (u) = det g T (w + 1) = det g <p{u + 1) is known as the quantum determi- 
nant of the quantum monodromy matrix. These formulas become more simply looking 
in the normalization (5.5): 

N 



W(u) = J2(-±) a T a (u)e 



W-\u) = J2 e ~ s9uJ s( u + !)■ 



(5.59) 



a=0 



s=0 



It would be interesting to clarify the role of the classical Lax operator (3.32) in the 
quantum spin chain. 



5.5 Undressing Backlund flow and Baxter Q-operators 

The construction of Q-operators originally introduced by Baxter [39] in his solution to 
the 8- vertex model is now generalized and developed in the literature in various directions 
(see e.g. [331 EH SI])- In this section, we explain how the Baxter Q-operators emerge 
from the master T-operator using the approach of [S]. 

5.5.1 Preliminary remarks 

We will construct a set of Q-operators which can be identified with the Baxter Q- 
operators in the sense that: 

1. The Q-operators are building blocks for the T-operators meaning that the latter 
are written in terms of the former in the Wronskian form. The Wronskian formulas 
resolve functional relations for the T- and Q-operators known as TQ-relations. 

2. The Q-operators satisfy bilinear functional relations (the QQ-relations). With the 
additional input of analytical properties of the Q-operators (which are polynomials 
in u in our case), they imply the system of Bethe equations for their roots. 

9 It is known that one can generate the T-operators for infinite-dimensional representations just by 
changing the expansion point of the non-commutative generating series 36, 37 . It will be very interesting 
to see how this phenomenon shows up for the master T-operator and to realize formulas for T-functions 
in as operators. 
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Our construction of the Q-operators is directly related to the undressing chain of Backhand 



transformations from section 4.2 In our approach, the nested Bethe ansatz is represented 



by this undressing chain, where the size of the matrices involved in the Wronskian-like 



solution is successively reduced, as in equation (4.36). As it was shown in section 4.2, 
the undressing procedure basically consists in picking the singular terms of the function 
T(u,t + [z^ 1 ]) (in this section it is an operator) at the poles at z = Pi,P2, ■ ■ ■ ,Pn- Since 
the poles are in general of a high order, there are different versions of this procedure 
according to the possibility to consider different terms in the singular part of the Laurent 
expansion. As it was argued in section 4.1.5, they are essentially equivalent and differ 
by certain normalization factors only. Below in this section we consider in some detail 
two distinguished choices: picking the coefficients in front of the highest and the lowest 
poles. 



is a 



As it was already mentioned, the master T-operator introduced in section 5.3.3 
polynomial r-function in the variable u while the r-functions from section 4 are poly- 
nomials multiplied by the overall exponential factor rij =1 p", so that the definition of 
the master T-operator differs from the r-function ( |4.13[ ) by the factor (detg)~ u . More 
precisely, let \u) be a basis of eigenvectors of the master T-operator and let (uj\ be the 
dual basis such that (oo\u') = 5^')- Then 

(u>\ T(u, t) \u) = (det g)- u r u {t- (5.60) 

where the r-function in the r.h.s. depends on the state According to this, the 

undressing procedure should be slightly modified in order to preserve polynomiality of 
the r-functions at each level. 

In this section we will work in the basis e a where the twist matrix g is diagonal: 
ge a = p a e a , a = 1, 2, . . . , iV: 

( 0,0,— ,0 ,1,0,— ,0). 

i-1 N-i 

A basis of eigenvectors of g® L — g (g) g <g) . . . (g) g is then given by 

V v ' 

L 

|aia 2 • • • a L ) — e ai <g) e Q2 <g) . . . <g) e aL 

N 

with eigenvalues p ai p Q2 ■ ■ -Pa L = J^pfS where Mj is the number of indices in the 

i=i 

vector \a1a2 . . . (Xl) equal to i. These numbers can be regarded as eigenvalues of the 
operator 



10 



Mi = (I® (fc_1) ® rru <8> I®( L " fc )) , where = diag(0, 0, ■ ■ ■ , 0, 1, 0, • ■ ■ , 0) = % 

(5.61) 



fc - 1 i-1 N-i 



10 Note that is a Cartan element en of gl(N) in the fundamental representation and Mj is its co- 
product A( i_1 ) (en). In addition, g can be written as g = exp(^^ =1 logpieu) — YiiLiPi" with respect 
to the above eigenbasis. 
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on the states \a\a 2 ■ ■ ■ oil), so one can write g® L = p^p^ 2 ■ ■ .p^- N ■ It directly follows 

N L 

from the definition that ^ M; = L I 0i and thus ^ M, = L. 

t=i i=i 
In pi [8], it is shown that the T-operators can be explicitly written through dia- 
grammatic expressions, and T(u, t) is then a combination of permutations and diagonal 
operators. As a consequence, the operators M, commute with all the T-operators. In 
fact this follows from the GL(iV)-invariance of the .R-matrix and means that the eigen- 
states \u) of the T-operators can be classified according to the set of quantum numbers 
Mi = M t (\u)) = H Mi \u) and 

N 

M9 |w> = [[Pi 

i=i 

For what follows it is important to note that Mi(\u)) is the multiplicity of zero at z = pi 
of the eigenvalue of the operator (zl — g)® L on the state \u). More precisely, the following 
operator identity holds: 

{zl- gf L = {z- Pl ) M ^z -p 2 ) M >...( z - PN ) M ». (5.62) 



5.5.2 The undressing procedure 



For our purpose in this section the most convenient version of the undressing procedure 
is picking the coefficient in front of the most singular term of the Laurent expansion 
of T(u, t + [z -1 ]) around its poles. As is shown in section 4.2, this version differs from 
the one based on taking residues by a normalization factor (which is an operator in the 
present section). 

Applying the undressing procedure based on the highest poles, we note that there 
are N different ways to define the T-operator at the first level of nesting by choosing 
i G {I,..., iV} and picking the coefficient in front of the highest pole at pf 



T®(u,t) 



■l) i - 1 det^pr 1 e -^ t 'P0 res 2 



-n [(z-p t ) M ^T(u+l,t + [z- 1 



(5.63) 



where Mj is the operator defined in (5.61). This "operator multiplicity" reflects the fact 
that the order of the highest pole depends on the sector of the Hilbert space where the 
diagonalization is performed and is equal to Mj + 1. This expression (5.63) is nothing 
but the equation (4.42) with an extra factor 



(0, 



u 



t)\u) 




which ensures the polynomiality of T®(u, t). Here we introduced a notation 
the state dependence of the function. 

As the simplest example, consider the case L 



T (L=1) (it, t + [z- 1 ]) = ((« - £i)I + kt k9 k + 



1. From (|5.29|) we find: 
9 



k>l 



Z 



(5.64) 



to denote 



e £;>i*itiV 

g J det(I - z- x g) ' 



(5.65) 
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As z — > Pi, the r.h.s. has a double pole in the subspace spanned by the basis vector 
and a simple pole when restricted to the subspace spanned by e a with a 7^ i. In the 
former subspace the eigenvalue of (z — Pi) Mi is z — pi while in the latter it is 1. The 



explicit calculation of (5.63) gives: 



N 



Pj 



^2 (" • 1 X! kl>i 'i l ' i>,< " 



k>l 



Pi-Pj 



l<a<7V 



p a e 



e(t,Pa) 



(5.66) 

where are the standard generators of the algebra gl(N). 

At the next level of the Backlund nesting chain, one should choose j G {1, . . . , iV}\{z} 
and define: 

TM(u,t) = (-l)^\det gfipip^e-^'^-^'^ipj - p { ) 

res 2i=pi [(zi —pi) Mi (zj —pj) M 'T (u + 2,t + fa 1 ] + [zj 1 ])} . 



(5.67) 



x 



z 3=Pj 



The sign factor is (— — where £y = 1 if i < j and —1 if z > j. 

With this sign convention (u,t) is symmetric in the indices i,j, or, equivalently, the 
transformations T — > T"W — > T™ and T — > T"^ — >■ Tw») commute for any 2, j, namely, 



y(y') — yw. The expression (5.67) is nothing but the n = 2 case of (4.43) with an extra 
factor 



N 



(w\TM(u,t)\u) 



n pri^t) 



(5.68) 



This factor is necessary to make r™(w, t) a polynomial of the variable u. 

In order to proceed to higher levels, one should fix a subset / = {21,22, • ■ ■ ,in} C 
{1,2,..., iV} and denote by J its complement: J = {1, . . . , iV} \ {ii, i 2 , . . . , i n }- The 
"nesting level" is defined to be the number n of indices in the set {11,12, ■ ■ ■ ,i n }- in 
complete analogy with the cases n — 1,2 considered above, we define T-operators 
y(ni2-i ra )( M; t) on all intermediate levels by the formula (cf. ( 4.43[ )): 



^^i„-l-Card{fc|i fe <i„,l<fe<n-l} ^ e ^.^ j J^Jpr 1 j e ~5(t,Pi n ) 

\a=l / 

xres 2!n=Pin [(^-pj^^-^ )(t* + l,t + [«- 1 ])], (5.69) 



which leads to 



T (n...*„) M) = (_i)^A n ( Pil , . . . ,pj(det <?)" n^" e " C(t,Pl ° 



va=l 



x res 



Kk<n 



n< 

ig/ 



Zi "Pi 



T(u + 71,1 + ^21^}) 



igj 



(5.70) 
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where A n is the Vandermonde determinant (4.37) and d n is defined in section 4.2. The 



T-operator at level n = is the master T-operator defined in (5.26). According to 



this definition, the correspondence of eigenvalues of the T-operators at level n and the 



r-functions at level n given by (4.44) is as follows: 



J6J 



(5.71) 



l<"> 



The r-function in the r.h.s. is obtained by the undressing procedure from r u (t; \uS)) 
according to the rules presented in section 4. Note that the eigenstate \ou) is common for 
all the T-operators. One can also introduce 'Thicker coordinates" for the T-operators 
on intermediate levels by the expansion 



T^ li2 - in \u,t) = ^s A (t)T (l1 



l2-.«n),A ( 



U 



(5.72) 



whose form is identical to (5.9). The diagrams A correspond here to representations of 



the subalgebra gl(N — n) C gl(N). 

It is clear that the bilinear relations of the form ( |4.39[ ), ( |4.40[ ) (cf. ( |5.19[ ), <\5.20h ) hold 
for the T-operators T^ x - ln \u^t) at any level and relate T^ 1 -^ 1 ^ (u, t) to T^ %1 - ink \u, t) 
and T^ 1 - in \u, t): 



Eijp^T^-^iu, i)T^- ink \u + 1, t) + e jk pT l T^--- i ^ k \u,t)T^--- ini \u + l 1 t) 

+ e ki p- l T^-- A - ki \u,t)T^--- i ^\u + l,t) = 0, 

p j 1 T( ix - i " i \u,t)T^- i ^(u + l,t) -pr 1 T( il - i »J)(M,t)T( <1 "- i » < )(u + l,t) 

= e iJ -r( il - i " i ^(M,t)T( <1 - i ")(M + l,t), 
where z, j, A; e{l, 2,..., JV}\{i 

We thus have the undressing chain that terminates at the level N: 

T(u,t) -> T (n) (M) T (ili2) (u,t) T^ X2 '" N \u, t) -> 0. 



(5.73) 



(5.74) 



(5.75) 



It can be shown that the last (the completely undressed) member of this chain is an 
invertible operator which does not depend on u, t: T^ 12, ^ N \u, t) = T^ 12 - N \ An advantage 
of the undressing at the highest poles is that this operator has a simple explicit form: 



A? 



rp{V2...N) 



det f 4- 

l<i,j<N 1 1=1 



(det g) N g® L 



N 



(5.76) 



Kl<k<N 



which follows from the comparison of the large u asymptotics of (4.13) and (5.26) 



Namely, from (4.13) and (4.9) it follows that 

N 



rjt) 



(dets)V^ (fl (a jMj e^) ^det^ (pjT"") + O (t*" 1 ) J (5.77) 
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while for T(u, t) we have 



T(u, t) = u L I® L eZti^Pu) + o (u 1 - 1 ) . 



(5.78) 



Then (5.76) follows from (5.77), (5.78) , (5.71) and the expression for the eigenvalues of 



the T-operator on the last level of nesting (follows from (4.44) and (5.71) for n = N): 



N 



H T( 12 -^ \u) = H(M j \a jMj 

3=1 



|w> 



Set 



Q j {u,t)^(T^--- N Y 1 T^--- N \u,t). 



In the notation of section 4 (see eq. (4.44)), we can write 

(u\ Qj(u,t) \u) = p] 



MA a 



(5.79) 



(5.80) 



where the r.h.s. should be evaluated for the r-function corresponding to the state \oj). 
Similarly, one can define 



Qj(u, t) = (T^ 2 -^)- 1 T^\u, t), J = {1, . . . , N} \ I . 



(5.81) 



In this general notation, the previously defined operator is Qj(u,t) = Qyy(u,t). From 
equation (4.44) we have, for J = {ji, j 2 , . . . , jN-n}'- 



In particular, at J = {1, 2, ... , iV} this yields 



(5.82) 



N 

T(u,t) = g® L (JJM* 

i=l 



det (p% +l - l Q k (u- I + l,t)) 



det {$>l 

Kk,KN K 



(5.83) 



where we have used (5.76). These Wronskian-like determinants can also be obtained by 



taking residues of (5.23) at z ik — >• p ik . The Jacobi identity for the determinant (5.82) 
gives the so-called QQ-relations, from which the Bethe equations follow. 

In fact the operators Qj(u, t) introduced in the way explained above are the (specially 
normalized) T-operators on the intermediate levels of the undressing procedure. The 
Baxter Q-operators Qj{u) are then defined as their restrictions to zero values of t: 



Qj(u) = Qj(u,t = 0). 



(5.84) 



We see from the above formulae that Qj(u) is expressed as a residue of the master T- 



operator. Moreover, if we restrict (5.80) to t = 0, and expand it similarly to (4.10) then 



we see that the coefficients a,i m in (4.9) can themselves be expressed from the residues of 
the master T-operator. 
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Comparing (5.11) to the equation ( 4.21[ ) (restricted to t = 0) and using the definition 
(5.72), we conclude that for the symmetric representation A = (s) = (s, 0,0, ...) the 



following formula holds: 



T (h...i n ),X {u) =T (12...N) ^ fX l+ l-l, 

KKN—n 



l Q h (u-l + Xi + l)j . (5.85) 
The Jacobi identity for this determinant is known as the TQ-relation. For all other 



representations this result can be plugged into an analogue of the CBR-formula (1.3): 



A'-l 



T («...in),A( u ) = f TT T^-^iu-k)) det T {h - in) ^- i+j \u-] + l), 

' i,J=l,.--,Ai 



(5.86) 



to get that the determinant formula (5.85) actually holds for arbitrary Young diagrams 



with at most N — n rows (if it has more than N — n rows, then the T-operator van- 
ishes identically). In the particular case 7 = equation (5.85) yields the Wronskian 



determinant representation of the T-operator at fixed representation: 



N 



T x (u) 



det [p N k +l - l+Xl Q k {u 

Kk,KN v 



i+Xi + r 



det (p£ 

l<k,l<N K 



N-l\ 



(5.87) 



At L = the pre-factors and the operators Qk(u) become simple (equal to l/pj) and 



this formula coincides with the determinant representation (3.21) for characters. By 



construction, all the Q-operators are polynomials in u and the determinant formula (5.87) 
at A = must give the fixed polynomial (p{u) = [7 



This requirement implies 
the system of Bethe equations for roots of the Q-operators (or rather their eigenvalues) . 

The undressing procedure is schematically illustrated in Fig. [2] for the case N = 



2. The green double arrow corresponds to the Backlund transformation (5.91), which 
"undresses" the T-operators by gradually decreasing the size of the matrices at each 
site of the lattice (which is the maximal possible height for the Young diagrams giving 
rise to non-zero T-functions). The blue arrow shows how the original T-operators are 
reconstructed by "dressing" with the Q-operators. 

5.5.3 On other versions of the undressing procedure 

The undressing procedure and the Q-operators can be defined in a number of equivalent 
ways. Here we mention some of them. 



Lowest poles. Instead of picking the highest order pole, one might consider another 
version of the undressing procedure, based on taking residues at other poles (not of the 



highest order), as in section 4.2. For example, at the first level (cf. (4.30)), we may define: 
T«M) oc res, =Pi [((^P^e^T^ + l,t + [z" 1 ])] . (5.88) 
This definition does not use the operators M i; at the price of less transparent normaliza- 



tion factors. In particular, the definition (5.88) with the unit pre-factor leads to a more 
complicated operator yl 12 - -^). 
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^(1,2). 




Figure 2: Backlund flow for the GL{2) spin chain. 



Construction of |8j. At last, let us compare our present construction with the one 



introduced in [8]. On the first level, noticing that (z — p t 



we can 



rewrite (5.63) as 



T«(w,t) = det^rV^IVes^ [(zl 



where N, = (-1)* 1 ] [ (p t - pj 



L T{u + l,t + [z~ 1 })] , (5.89) 

(5.90) 



-M,- 



The normalization factor N, is independent of u and t, and commutes with all the T- 
operators. For this reason it was disregarded in [8]. 



On higher levels (5.70) can as well be written as: 



n 

r (H-*n)( w ,t) = (detsr(n^ e ~ f(t ' Pla) ) Ntl - 



a=l 



x res 



Kk<n 



.7=1 



a=l 



(5.91) 



where N ix An is an operatorial normalization factor which we don't specify here. Let us 



now show that (5.91) coincides with equation (32) from [HJ. 
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1. In [S], the set of times t was not introduced. But as it can be seen from (5.31), the 

r> .. . . .r l-1. . i . . . i . .. r> / T \ /i./tt 



transformation t t- >■ t + [z l ] is equivalent to the insertion of w(z x ) = (det(I 
z~ 1 g))~ 1 to the right of the co-derivatives. 

2. The limit Zj — > l/a;^ used in [5] has the same effect as taking residue at the point 
Zj = pj if we change notations z — > 1/z, pj — > Xj, and notice that res z=p f(z) = 
lim z ^ p (z — p)f(z) for a function / with a simple pole. 

3. In [8], depending on the representation, an overall factor is necessary to convert 
characters of g = diag (Pi\i <i<N ) G GL(N) into characters of gj = diag fpilj 6 jl £ 
GL(| J|). In our current notations, this amounts to changing F g (t) to F gj {t). One 
can easily see that this is exactly what the factor e~^<*=i€( t >Pia) does. Indeed, 
e-^^ a ) = Fgj{ t)/F g (t). 



This comparison shows that the undressing procedures of [8J and the one of the present 
paper are essentially equivalent. But in the setting of our present paper this undress- 
ing procedure has a transparent connection to the classical integrable hierarchy and its 
polynomial r- functions. 

Our construction clarifies the combinatorial structure of the set of the 2 N Q-operators 
on the Hasse diagram found in [32] • Moreover, it naturally takes into account the fact 
that all these Q-operators, as well as all the T-operators, are polynomials in the spectral 
parameter u. As a consequence, the Bethe equations are automatically contained in this 
analysis, as it was explicitly demonstrated in [8J. 



5.6 Dynamics of zeros of the master T-operator 



As we have seen, for a GL(iV)-invariant generalized spin chain on L lattice sites the master 
T-operator T(u,t) is essentially a polynomial in u of degree L. It can be represented in 
the form 



T(u,t) oc Y\_( u ~~ u j(^)) 



(5.92) 



where the roots Uj depend on all other times £j (we have omitted an irrelevant exponential 
function in front of the product). In some special case; 



ii 



one or more roots Uj can tend 



to infinity and so the degree of the polynomial can be less than L. Note that we treat 
equation (5.92) as an operator relation and the roots Uj(t) are (commuting) operators. 



Another possibility is to treat (5.92 ) as a relation for eigenvalues of the master T-operator, 



then the roots of each eigenvalue have their own dynamics in the times tj. This dynamics 
is known [19] to be given by the rational Ruijsenaars-Schneider model. A solution of the 
spectral problem for the master T-operator in terms of the integrable dynamics of its 
zeros and a connection with Bethe equations will be given elsewhere. As is well known, 
the Ruijsenaars-Schneider model contains the Calogero-Moser system of particles as a 
limiting case. In this connection let us note that a similar relation between the quantum 
Gaudin model (which can be regarded as a degeneration of quantum spin chains or vertex 



ii 



For example, when the boundary twist parameters tend to 1 
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models with rational .R-matrices) and classical Calogero-Moser system was found in 
from a different reasoning. 



An important remark is in order. As is seen from (5.3 ), the inhomogeneity parameters 
£i are coordinates of the Ruijsenaars-Schneider particles at U = 0: & = itj(0). Clearly, 
they provide only a half of the initial data necessary to fix the dynamics. The other 
half is initial values of momenta which are implicitly determined by the above mentioned 



conditions on the r-function (5.9). Since the quantum transfer matrix for the spin chain 
is of the size N L xN L , it has not more than N L different eigenvalues. This means that the 
system of equations that fix the initial momenta has only a finite number of solutions, 
corresponding to common eigenstates of the transfer matrices. 



6 Concluding remarks 

Our main result in this paper is the identification of the generating function of com- 
muting integrals of motion in generalized quantum spin chains (called here the master 
T-operator) with the r-function for integrable hierarchies of classical soliton equations. 
Among other things, this result means that quantum integrable models have a hidden 
free fermionic (Grassmann) structure in the auxiliary ("horizontal") space. It would be 
very interesting to relate it to the hidden Grassmann structure uncovered in H5] in 
the quantum ( "vertical" ) space of quantum integrable models. Another context where r- 
functions of classical integrable hierarchies enter the theory of quantum integrable models 
and associated 2D lattice models of statistical mechanics is calculation of scalar products 
and partition functions with domain wall boundary conditions [IHJ H3 08] . 

Many things remain to be done in order to further clarify and generalize this classical- 
quantum correspondence and, possibly, to use it for inventing new methods to solve both 
kinds of problems. Here we list some obvious directions of further development: 

- To construct the Bethe vectors within this formalism (which does not use the 
standard algebraic Bethe ansatz) and to clarify their role in the classical MKP 
hierarchy and its polynomial reduction according to Krichever. 

- To extend the results to quantum integrable models with trigonometric i?-matrices. 
This is more or less straightforward and clear how to do. Whether the approach 
based on the co-derivative is applicable is not quite clear at the moment, however. 
Instead of rational solutions of the MKP hierarchy one should consider solutions 
with trigonometric dependence on the spectral parameter. The Krichever's con- 
struction is directly extendable to this case. Zeros of the master T-operator move 
as particles of the trigonometric Ruijsenaars-Schneider model. 

- To extend the results to quantum integrable models with elliptic (Belavin's) R- 
matrices. This is the most general and most interesting case. Zeros of the master 
T-operator move as particles of the elliptic Ruijsenaars-Schneider model. The rel- 
evant class of MKP solutions is a special subclass of general elliptic (two-periodic) 
solutions. One of the interesting problems here is to characterize this subclass 
in algebro-geometric terms through the corresponding algebraic curves. Elliptic 
solutions in general position correspond to smooth curves which are coverings of 
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elliptic curves. Presumably, the solutions relevant for quantum spin chains are not 
of general position and the corresponding algebraic curves are singular. 

- To extend the results to supersymmetric GL(N |M)-invariant integrable spin chains. 
The general definition of the master T-operator and the Hirota equation remain 
the same. The realization through the co-derivative is known from j8], and can 
certainly be related to Krichever's rational r-function. 

- The fact that the Baxter Q-operators and T-operators on intermediate levels of the 
nested Bethe ansatz are contained in the master T-operator seems to be general 
and not restricted to models with rational i?-matrices. One can try to find a 
general procedure to obtain these objects from the master T-operator applicable in 
all the cases listed above. The example of rational solutions suggests that such a 
procedure should be formulated in terms of Baker- Akhiezer functions on singular 
algebraic curves. 

- A development of this work might also contribute to representation theory of quan- 
tum affine algebras and, in particular, to the theory of g-characters introduced by 
Frenkel and Reshetikhin [35] as a spectral parameter dependent version of usual 
characters. 

- We also expect that our result can be extended to quantum models associated 
with more general Lie algebras. This is in accordance with the fact that the KP 
hierarchy has analogs for Lie algebras other than A^-i [16J. 

- Our classical interpretation of the quantum integrability should be also applicable 
to the 2-dimensional integrable quantum sigma models, where the operatorial for- 
malism is almost unexplored. The most interesting case of such applications could 
be the construction of an operatorial form of the AdS/CFT Y-system (and the 
associated T-system) [37J 113 EH EE] and clarification of its relation to the classical 
integrable hierarchies. 

A more ambitious goal is to suggest an alternative approach to diagonalization of 
quantum transfer matrices which would not use the Bethe ansatz at any step. This does 
not seem hopeless because the quantum problem appears to be reducible to a purely 
classical integrable dynamics of Calogero-like systems of particles. 
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Appendix A: The generalized Wick's theorem 



Let Vi = J2j v ij' l Pj be linear combinations of ^-'s only and w* = Yuj w lj' t / } j be linear 
combinations of ?/>*'s only. The following version of the Wick's theorem holds: 



(n 


V\...V m 


«4 


...w\G 


n) 


(n 


G 


n) 



det 



(n 




KG 


n) 


(n 


G 


n) 



(Al) 



where G is any group-like element. This statement can be proved by induction using the 



bilinear identity (2.22). Indeed, suppose (2.23) is valid for some m > 1 (it is trivially 



valid at m — 1). Set 

(U\ = (n\w* 1 , (U'\ = (n\v 1 v 2 ...v m+l w* rn+1 w* m ...w* 2 i \V) = \V) = \n) . 

Plugging this in the bilinear identity, we see that its r.h.s. vanishes while the rest yields, 
after some transformations: 

(n\ G \n) (n\ V\ . . . t> m +iu>„ +1 . . . w*G \n) 

m+l 

= (n\ Vjw\G \n) (n| v x . . . fa . . . v m+ iw* m+l . . . w* 2 G \n) 

3=1 

(here fij means that Vj is omitted). But the r.h.s. is the expansion of the determinant 
for m + 1 in the first column. 



To show that the two determinant formulas (2.24) and (2.25) are equivalent, let us 



take 



(E/| = <n + l|Vz<, (U'\ = (n\, |V> = |n>, \V') = \n + l) 



(A2) 



in the bilinear identity (2.22). Using the relations (2.6)-(2.9), we see that only one term 
in the r.h.s. remains: 

(n + l\ipiW*ip k G \n) (n\ ?p* k G \n + 1} = (n + 1| ipiW*G \n + 1) (n| G \n) 

k 

In the l.h.s. we plug w* = ^ Kj^j an d move ipk in the first factor to the left position. 
Using the anti-commutation relation for the fermion operators, and transforming the sum 
J2k Kk^Pt back to w* in the second factor, we arrive at the 3-term identity 

(n\ G \n) (n + 1| ipiW*G \n + 1} = (n + 1| G \n + 1) (n| ipiW*G \n) 

+ {n+l\ipiG\n)(n\w*G\n + l). (A3) 
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Let us divide it by (n| G |n) (n + 1| G |n + 1), take a sum over the values of n equal to 
n — j, n — j + 1, . . . , n — 1, and then put I = n — j. Then we obtain the relation 
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n) _ {n-j 


w*G 


n-j+1) 
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n-j +1) 
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fc=i 



(n— A; +1 V'n- 




n — k) (n — k 


w*G 


n-k + 1) 


(n-Jfc|G 


n-k) (n-k+l\G 


n-k + 1) 



(A4) 



which shows that the columns of the matrix in (2.24) are the linear combinations of the 



columns of the matrix in (2.25) and their determinants are thus equal. 



Appendix B: Proof of the quantum Jacobi-Trudi for- 



mulas for the Pliicker coordinates (3.18) 



In this appendix we prove the "quantum Jacobi-Trudi" formulas (3.18) for the coefficients 



c x (n) (3.16) 



Let us begin with the case when A = (ac\(3) is a hook. We have: 
c {a \p)(n) = (-1)^ +1 (n\ ^Vn-ZJ-iG \n) 

= (n-p-M r n - P -i ■ ■ ■ r n -ir n+a i>n-p-iG \ n ) 

= (n-{3-l\r n -p..-r n -ir n+a G\n). 



Applying the Wick's theorem in the form (2.25), we get 

, , (n-j\iJt-i +n G\n-j + l) 

det — -, . \ , „ , ; — 77 — , 

j,j=i,...,/3+i {n-j + l\G\n-j+l) 

where Ai = a + 1, A2 = A3 = . . . = A^+i = 1. The numerator in the r.h.s. can be easily 
identified with C(Xi-i+j-i\o) {n ~ J + 1)> so we finally obtain 



{n-P-1 








n) 


(n 


G 


n) 



c {a \p){n) = r n (0) . . det — 



(Bl) 



TT c (n - k) det c Xl - i+j (n -j + 1), 

t=i / ' 



where c s (n) := C( s _i|o)(n) are the expansion coefficients for one-row diagrams. 

In a similar way, one can prove that the determinant formulas of the Jacobi-Trudi 



54 



type hold for arbitrary Young diagrams. Let us rearrange (3.16) as follows: 



d(X) d(X) 

c x (n) = („| Hr n+ai H^-iG \n) 

i=i i=i 

d(A) d(A) 

i=i i=i 

ft + l d(A) d(A) 

= (-i)HA)^(A) + ^^ ( n - & - ii nvc-*n^-iii^ G i"> • ( B2 ) 

fe=l j=l i=l 



By using the relations ip^m — 1 — VvVCi an d (w— ft — 1| V'm = for m G {n— ft — , 
we obtain 



S+l d(X) 

c x (n) = (n - ft -1| II^UlI^GIn). (B3) 

Taking into account the following relations for the Frobenius notation: {i — Ajjfl^ = 

{-aii}t=i, - ^i}fid(A)+i = {1, 2, . ..,/?! + 1} \ {ft + we can rewrite this as 

follows: 

ft+i 

c A (n) = (n-ft-l| n^+A^N- (B4) 

i=i 

Applying the Wick's theorem in the same way as for the A = (a\/3) case, we finally obtain 



the desired result (3.18). 



Appendix C: Proof of the quantum Giambelli formula 

The quantum Giambelli formula was proposed in [32] in relation to the analytic Bethe 
ansatz for finite dimensional U q (B^') modules. It was also mentioned there that this 
formula follows from Sylvester's theorem without detailed calculations. To make the 



text self-contained, we will give a proof of (5.7) based on the CBR formula (5.6). Here 



we follow [201 P a g e 88, example 22], where the classical case (without spectral param- 
eter) is discussed. In this section, we will use Frobenius notation for a Young diagram 
mentioned in section 2.2. We introduce the following notation: t^j = Tj_i(u — j + 1), 
n = \[, d = d(X), 4>{u — k + 1) = Tq{u — k + 1) = t kjk . We also introduce tuples 
I = (—oil, . . . , —a^ 1,2, ... ,n), J = (1,2, ... ,n). We will use the following relations for 
the Frobenius notations: {i-Xi}f =1 = {—an}f =1 , 0"-A*}f =d+1 = {1, 2, . . . , n}\{ft + l}f =1 . 



Then (5.6) can be written as 

T\u) = (f[t k>k ) det {U- Xi)j ) = f f\t kk ) det fa). (CI) 
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Next we introduce the following matrices: 

Mi = (*ij)(ij)6iixJ) M 2 = (ti,j)(i,j)eJxJ-> U = (5j )( g i+ i)(ij) 6 j x j i; (C2) 

A = ^ 17 M 2 J ' B = { U M 2 ) ' (C3) 

where 7 X = (— ckj, . . . , — a^), Ji = (1, 2, . . . , <7) , and Irf X ^ is a <i x <i unit matrix and (0)dxn 
is a d x n zero matrix. Then we find 

det = {-VpU ^ +d2 detB. (C4) 

Let us replace j-th row of A with i-th row of B and write the determinant of this matrix 
as Cij. We will consider a matrix C = (cij)i<ij< d+n . Taking note on the relations 
Cij = Sij det A for d+1 < i,j < d + n and = for d+l<i<d + n, 1 < j < d, we 
obtain 

det (q,)= det (cj ,)(det A) n . (C5) 

l<i,j<d+n ,J l<i,j<d 



We will use det A = YYj=i Expanding the j-th. row of the determinant Cjj, we obtain 

n+d 

f * l<p,q<n-{-d, 

k=l p^j.i^k 
n+d 

= KkiA- 1 )^ det A = {BA-\j det A (C7) 
fc=i 

Then we obtain an identity 

det (a^ = det5(det A) d ~ l . (C8) 

l<i,j<d ' 



We can calculate (CI ) for the hook diagram A = (a + 1, b): 

/b+l \ _1 6+1 



l<p<b,l<q<b+l, 

\k=2 / k=l q?k 



For 1 < i, j < d, we can calculate explicitly 

n 

dj = E(-l) ft+/ W .det (W) (CIO) 

fc=l p£pj+l,q£h 

Pj + l n 

.^(-i)^^ n ^ ten) 

fe=l g^fc p=/3j+2 

n 

= (-1) A+ Xa(«)II*w». ( C12 ) 

p=2 
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where we used (C9). Based on this identity, we obtain 



p=2 



Hence, from Q, (|UT}, Q, flCT3] ) we get 

T X {u) = (ii,i)~ d+1 det (Tq,. «. (m)). 



l<M<d 



This is (5.7). 



(C13) 



(C14) 



Appendix D: The co-derivative 

In this appendix we summarize formulas on the co-derivative introduced in [9]. 

The (left) co-derivative is the operator that acts on any function / of a group element 
g E GL(N). Symbolically, it is defined as follows: 



D®f{g) 



d 



® f(^ e g) 



(Dl) 



</>=o 



Here 



12 



N 



■ e 



XI 0fl e a/3 is a general element of the Lie algebra gl(N) parametrized 



a,8=l 



by (fig, with e a p being the matrix with only one nonzero (equal to 1) matrix element at 
the place (a/3): (e a p)% = S afl 8p u . To put this more precisely, let us start from scalar 
functions. In this case, we do not need the symbol £g> and get: 



Dtrg r 



ng' 



D det g = det g ■ I, 



(D2) 
(D3) 



where I is N x N unit matrix. If / is a iV x iV matrix valued function, the result of (Dl ) 
belongs to iV 2 x iV 2 matrix. In particular, if f(g) = g, we have in components: 







- (e^g) : 

i V vjj. 



4>=o 



In invariant form, this formula can be written as 

D ® g = V(l ® g) 



(D4) 



where V is the permutation operator defined as V(A ® B) = (B ® A)V for any N x N 
matrices A, B. Explicitly, we have V = J2ab=i e ab ® e ba- 

In the main text, we considered actions of the co-derivative on functions in (C N )® L . 
In this case, it is convenient to define the co-derivative in the following way. Let us 



12 In this paper, we assign the left matrix suffix to the superscript and the right suffix to the subscript. 
For example, element Ajj of a matrix A — (A^j) is written as Aj. 
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consider an embedding of any matrix A on into (C N )® L , and introduce the following 
notation: 



Ai = I ® ■ ■ ■ I ®A ® I <g> • • • ® I . 

i-l 



(D5) 



The action of Di on any function of g £ GL(N) in (C 



L-i 

N\®L ig defined as 



N 



d 



(D6) 



a,b=l 



where e^) = ^ £g> e a & (g) I®^ L ^. In components, (D6) reads 



(A/^))S:S = Es|r(/( e ^)) 



Cn...«i_l7«i + 1...Q!L 



/9i /8i 



4>=o 



For example for /(g) = gr^, we obtain a generalization of (D4): 



(D7) 



where Vij = ^2 ab e i& e &a- For i = J case, the permutation operator reduces to the 
second order Casimir, and in this case, it is Vu = ^2 a b e ab e bl = ^2ab e ^ = NI® L since 
we are considering the fundamental representation. Then we obtain = Ngi. In 



particular for L = 2 case, D x g 2 = V\ 2 g 2 stands for (D4), while D 2 gi = V 2 \g\ stands 



for (I <g> D)(g (g> I) = V(g (8> I). Now the Leibniz rule for any functions A(g),B(g) of 
g £ GL(N) can be written as 



I=JUK,JnK=® 



\{biA{g) 



(D8) 



where I = (i u i 2 , . . . , i n ), J = (i ai , i a2 , . . . , i ak ), K = i A , . . . , i Pn _ J («!<■••< 
ttfc! /?!<•••< /3 n -fc) are, as sets, subsets of {1, 2, ... , L}, which satisfy 7 = J U K, J D 
K = The ordered product on any indexed operator {v4j}j £ / is defined as ffig/Aj = 
Ajj • • • Ai n . The Leibniz rule can be used, for example, 



n-l 



n-l 



n-l 



A# = E^^ k ; = E^^K"*" 1 = E 7 ^. 



n—k 
3 ' 



(D9) 



fc=0 



fc=0 



fc=0 



where n £ Z> . Let us introduce commutative element acting in (C ) 



N\®L 



c f } e n 

J2i£i m,=n,m i eZ> ig-T 



(D10) 



where 7 C {1, 2, ... , L}, n £ Z> , Cf } = I 0L and = for n < 0. Note that this 

(Dll) 



commutes with any permutation indexed by the I. ( D9 ) can be written as 



58 



This can be generalized for example as 



DiDjfi = C%-%{b igj ){b i9k ) + C^Dibjgn 

= C { ;J } V lj9j V lk g k + C^ } V jk V ik9k (D12) 

Let S (n) (J) be a sub group of the permutation group over the set J whose elements have 
exactly n cycles: any element a of S^ n '(J) has a cycle decomposition a = T\ • • -r n ({r.i} 
are mutually disjoint cyclic permutations on J). Then we have the following formula 

L y 

i=l luJ={i,2...,L} iel jeJ 

Card(J) f. T ^(*- ff -'(lfc)) 

= e n«. e «- e p.n ^' «w- 

/uj={i,2...,l} ie/ n=l o-eS(™)(j) feeJ y 

(D13) 

Here V a is a matrix in the space (C Ar )® L , with the matrix elements 

L' CT bl-Ji U «J CT (1) U «LJ CT (L)- 

The action of cr e S^ n '(J) on i 6 {1, 2, . . . , L} \ J is formally defined as a(i) = i. The 
summations are taken over all the decompositions of {1, 2, ... , L} into two disjoint sets 
/ and J. The above identity can be proven by induction in L and using the Leibniz rule. 



Yang-Baxter equation and commutation of co-derivatives One can check that 
these co-derivatives obey the following commutation relation 



bibj - bjbi = (A - b^Vij = -7> 4i (A - 

where the last equality comes from DjVij = VijDi- 

This relation can be understood from the following Yang-Baxter equation 



(D14) 



j-Ui + Vij) \\{u k + D k ) = Y[( U r{k) + b r{k) ) (Uj - Ui + Vij) 
k=l k=l 

{% if k = j 
j if k = i and % < j . 

k otherwise 



(D15) 



(D16) 



For instance, if we write (D15) when L = 2, and v>2 = 0, then the coefficient of degree 
one in u\ gives exactly (D14). Multiplying (D15) by Vij from the left, we can rewrite 
(D15) in the notation of the tensor product used in the main text: 

L L 

(1 + ( Uj - Ui)Vij) + D) = 6§{u T[k) + D) (1 + { Uj - ujVij) . (D17) 



k=l 



k=l 
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The equation (D15) can either be proven "from scratch" by proving explicitly the 
property (D14) of co-derivatives, or it can be deduced (when the co-derivatives act on 
linear combinations of characters, which is sufficient for this paper) from the usual Yang- 
Baxter Identity. 



Appendix E: Proof of (5.39) 



The left hand side of (5.39) is a polynomial of u±, u 2 , ■ ■ ■ , Ul- The degree of this polyno- 
mial with respect to each ttj (1 < i < L) is 0, 1 or 2. Taking note on this fact, we will 



prove (5.39) using the induction with respect to the length of the chain L. The proof 



that ( |5.39 ) holds for L = Lq, consists in the following three steps : 



(i) First we will see that the terms of degree two in one of the variables Ui are zero 
under the hypothesis that (5.39) holds for a chain with length L = Lq — 1. 

in Uj is obtained by the 



This step is easy since the term of degree two 



substitutions V Lo (u 
\\ ^ (CM 



Ui + 1 + D) <g) Ujl ® (g)f° 



Ujl 



■■3+1 



[Ui + D) and 
'ui + 1 + D). We see that the op- 



i=j+l 



erator obtained from the action of these co-derivatives is a trivial operator on the j th 
space times an operator acting on a spin chain with Lq — 1 spins (corresponding to all 
sites except the site j). We can then recognize that the coefficient of u 2 - in the left-hand 



side of (5.39) at L = L is equal to the left-hand side of (5.39) at L = L — 1 (up to a 



relabelling of the spaces), which is zero by the recurrence hypothesis. 

(ii) The other terms have degree one in some variables , Ui 2 , . . . , Ui n and degree zero 
in the other variables Uj 1 ,Uj 2 , . . . } Uj L n (where {ji, . . . ,jL -n} denotes the complement 
of {ii, . . . , i n })- We can show that all the terms where 1 G {zi, . . . , i n } vanish, assuming 
that ( |5.39[ ) holds for a chain with length L = L — 1. 

To perform this step, we have to investigate the coefficient of degree one in m, in 



each term of the left- hand- side of (5.39). For instance, in the operator 



V 



L 



u+lM^- 1 )) (v Lo (u)w(z^) 



the term of degree one in u\ is given by 



l \ / l q 

(1 + D) <g> + 1 + b)w{z^) I uil ® + b)w(z 2 l ) 

i=2 J \ i=2 

L \ / L 

+ I uil (8) (g)(«i + 1 + b)w(z^) lb®0( Ui + b)w(z 2 r ) 

i=2 J \ i=2 



Ml(I + D) ® 



L \ / Lq 

Ui + l + D)w{z^) (g)( Wi + D)w(z 2 l ) 

i=2 / \ i=2 




(El) 



where the right-hand-side is obtained by the Leibniz rule. The same analysis 
is easily performed for the other terms ( V Lo (u + Vjw^ 1 )) ( t> Lo (u)w(z^ 1 )) and 
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V Lo {u + l)w(z{ 



V 



, (appearing in (5.39)) and we obtain that the co- 

D) on the left-hand-side of 



efficient of degree one in u\ is obtained by the action of (I 
the identity with L = L — 1 spins. Hence the recurrence hypothesis shows that they 
vanish. 

(iii) Finally, we will deduce the vanishing of the term of degree one in the 



variables u^ , Ui 2 , 



and degree zero in the variables Uj 



3i i u h i 



where 



{ji, . . . , jL -n\ denotes the complement of . . . ,i n }. We will prove this vanishing 
by recurrence over n (i.e. we will assume that it holds for n — 1, and we will deduce it 
for n). 

First, let us assume that the terms with degree one in Uk and degree zero in u^+i 
(i.e. the terms where k G {zi, . . . ,i n } and k + 1 G {j±, . . . , j£ _ n }) do vanish, and let 
us deduce that the term with degree zero in Uk and degree one in uu+i (i.e. the terms 
where k G {ji, . . . , jho-n} and k + 1 G . . . , i n }) also vanish. To this end, we can use 



a Yang-Baxter equality (D17) for each term in the identity (5.39). For instance for the 
first term (i.e. for Z2\T) Lo {u + l)w(z^ 1 )) (V Lo (u)w(z2 1 ))) we write 



z 2 (l+(«fc+i - u k ) V k ,k+i 




(ui + l + b)w(zl 



- ^ ( (gx^+i+DM^r 1 

.1=1 





Ui + D)w(z 2 v 



u„u)+b)w{z 2 x ) (l + (u k+ i - U k ) V k ,k+l) 



(E2) 



where 



k+1 if i = k 
k if i = k + 1 
i otherwise 



(E3) 



If we sum this Yang-Baxter identity for each term in (5.39) and we keep only the 
terms of degree one in the variables , Ui 2 , . . . , Ui n and degree zero in the variables 



U jl 5 U 32 5 



u 



3 Ln — n. ? 



then we get an identity for the sum of the following terms 



First, the terms which contain Vk t k+i- As the total degree in u^ is one and the total 
degree in Uk+i is zero, these terms do not have any occurrence of u^ or u^+i in the 
factors containing co-derivatives. One can easily see that (up to multiplication by 



Vk,k+i) these terms are exactly the terms of (5.39) with degree zero in the variables 



U jl ) U 32 1 



i u 3L -n an d u k, and degree one in the other variables. By the hypothesis 



of the recurrence over n, these terms vanish. 



Next, the terms of the first line in (E2) which do not contain Vk,k+i- If we sum 



them for the three terms of (5.39), we exactly obtain the terms of degree one in 

Uj, of 



the variables u ix , u i2 , 



( [5T39] ) (with k G 
terms do vanish. 



u in and degree zero in the variables Uj x , Uj 2 , . . . , Uj Ll 
,i n } and k + 1 G {jt, . . . ,j Lo -n})- By hypothesis these 
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Finally, the terms of the second line in (E2) which do not contain Vk k+i- If we sum 



them for the three terms of (5.39), we exactly obtain the coefficient of degree one 
in Uk+\ and degree zero in Uk (and with degree one, resp zero in the other variables 
i m resp j m ). These last terms vanish (because the sum of expressions of the form 



(E2) is zero, and because the other terms vanish). Here we used the fact that 
the coefficient of u\u° k+1 in a polynomial f(u a ^,u a (2), ■ ■ ■ , Ua(L)) coincides with the 
coefficient of u° k u\ +l in a polynomial f(ui, u 2 , . . . , ul), where f(u%, u 2 , ■ ■ ■ , u£) is 
any polynomial in Ui, u%, ■ ■ ■ , ul- 

This shows that if the terms with degree one in Uk and degree zero in Uk+i (i.e. the 
terms where k G {ii, . . . , i n } and k + 1 G {ji, . . . , jio-n}) do vanish, then the term with 
degree zero in w fc and degree one in w fe+1 (i.e. the terms where k G {ji, . . . , jL -n} and 
k + 1 G {zi, . . . , i n }) also vanish. 

The same result is easily obtained if replace k + 1 with k — 1, and by iterating this 
argument, we obtain that the term of degree one in any set of n variables u^, Ui 2 , . . . , Ui n 
(and degree zero in the other variables) do vanish. Indeed, the above argument allows to 



reduce these term where 1 G . . . ,i n }, treated above. This proves the equation (5.39) 
by recurrence. 
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